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ABSTRACT 

In  this  paper  we  consider  the  following  problem:  Given  a  set  Q  of  n  (possibly 
intersecting)  line  segments  in  the  plane,  preprocess  it  so  that,  given  a  query 
ray  p  emanating  from  a  point  p,  one  can  quickly  compute  the  intersection  point 
$(^,p)  of  p  with  a  segment  of  Q  that  lies  nearest  to  p.  We  present  an  algo- 
rithm that  preprocesses  Q,  in  time  0(n'^''^  log'^  n),  into  a  data  structure  of  size 
0(nQ(n)  log"*  n),  so  that  for  a  query  ray  p,  ^(Q,p)  can  be  computed  in  time 
0{\/na(Ti)log^  n),  where  cj  is  a  constant  <  4.33  and  a{n)  is  a  functional  inverse 
of  Ackermann's  function.  If  the  given  segments  are  non-intersecting,  the  stor- 
age goes  down  to  O(nlog^n)  and  the  query  time  becomes  0{y/n\og^  n).  The 
main  tool  that  we  use  is  spanning  trees  (on  the  set  of  segment  endpoints)  with 
low  stabbing  number,  i.e.  with  the  property  that  no  line  intersects  more  than 
0{\/n)  edges  of  the  tree.  Using  such  trees  we  also  obtain  ftister  algorithms  for 
several  other  problems,  including  implicit  point  location,  polygon  containment 
and  implicit  hidden  surface  removal. 


^Work  on  this  paper  has  been  supported  by  Office  of  Naval  Research  Grant  N00014- 
87-K-0129,  by  National  Science  Foundation  Grant  DCR-83-20085,  and  by  grants  from 
the  Digital  Equipment  Corporation,  and  the  IBM  Corporation. 

'A  preliminary  version  of  this  paper  appears  in  Proceedings  5"*  Annual  Symposium 
on  Computational  Geometry,  1989,  pp.  315-325. 
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1      Introduction 

In  the  last  few  years  many  efficient  randomized  algorithms,  based  on  the  random  sampling 
techniques  of  [CI]  or  on  the  related  e-net  theory  [HW],  have  been  developed  to  solve 
efficiently  a  variety  of  geometric  problems.  One  such  recent  development  is  due  to  Welzl 
([We];  see  also  [CW]),  who  showed  that,  for  a  given  set  S  of  n  points  in  the  plane,  there 
exists  a  spanning  tree  T  of  5,  such  that  no  line  intersects  more  than  0{y/n\ogn)  edges  of  T . 
Such  a  tree  T  is  called  a  spanning  tree  with  low  stabbing  number  (a  formal  definition  is  given 
in  Section  2).  Welzl  used  spanning  trees  with  low  stabbing  number  to  obtcdn  an  almost- 
optimal  algorithm  for  simplex  range  searching,  namely  given  a  set  S  of  n  points  in  the 
pleine,  preprocess  it  into  a  data  structure  of  linccir  size  so  that,  for  a  query  triangle  A,  one 
can  quickly  count  (or  more  generally  report)  all  points  of  S  lying  inside  A.  His  algorithm 
counts  (resp.  reports)  the  points  lying  inside  a  query  triangle  A  in  time  0(>/nlog^n) 
(resp.  0{y/n\og^  n  +  K),  where  K  is  the  number  of  points  inside  A).  Soon  after  this 
paper,  Edelsbrunner  et  al.  [EGH*]  used  these  trees  to  preprocess  a  given  set  £  of  n  Unes 
in  the  plane  into  a  data  structure  of  size  0{n  log  n)  so  that,  for  a  query  point  p,  the  face  of 
the  arrangement  AiC)  containing  p  can  be  computed  quickly.  The  main  challenge  in  both 
of  these  papers  was  to  use  only  roughly  linear  space  (i.e.  0(n  log  '^'n)  space),  because 
if  we  allow  quadratic  space,  then  a  query  can  be  easily  answered  in  O(logn)  time  [Ed], 
[EOS],  [EG]. 

In  this  paper  we  present  several  new  applications  of  spanning  trees  with  low  stabbing 
number.  The  algorithms  presented  in  this  paper  are  faster  than  the  previously  best  known 
algorithms  for  these  problems.  One  of  the  main  goals  of  this  paper  is  to  demonstrate  that 
such  a  spanning  tree  is  a  versatile  tool  that  can  appHed  to  obtain  efficient  algorithms  for 
a  large  class  of  problems,  much  beyond  the  simplex  range  searching  problem,  for  which 
they  were  originally  introduced.  We  also  show  that  by  combining  the  spanning  tree  data 
structure  with  the  recent  partitioning  algorithm  of  [Aga],  [Agb],  we  can  reduce  the  query 
time  if  we  allow  more  space.  Similar  tradeoffs  between  space  and  query  time  have  been 
obtained  earher  [EGH*],  [Age],  [Chd]. 

The  first  and  perhaps  the  most  interesting  appHcation  that  we  consider  is  ray-shooting 
in  arrangements  of  segments.  There  are  two  versions  of  this  problem,  one  for  segments 
that  axe  non-intersecting,  emd  one  for  an  arbitrary  collection  of  segments.  Formally,  these 
problems  can  be  stated  as  follows: 

(a)  Given  a  collection  ^  =  {ci,  . . .  ,  e„}  of  n  non-intersecting  line  segments  in  the 
plane,  preprocess  it  so  that,  given  a  query  ray  p  emanating  from  a  point  p  in 
direction  d,  we  can  quickly  compute  the  intersection  point  ^{Q,p)  of  p  with 
the  segments  of  Q  that  lies  nearest  to  p  (see  figure  la). 

(b)  Same  problem,  except  that  the  segments  in  Q  can  intersect  arbitrsLrily  (see 
figure  lb). 
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Figure  1:   Ray  shooting  in  an  arrangement  of  (a):  non-intersecting,  (b):  arbitrary  segments 


If  the  segments  in  Q  form  the  boundary  of  a  simply  connected  region,  then  the  algorithm 
of  Cha^elle  and  Guibas  [CGa]  preprocesses  Q  into  a  data  structure  of  linear  size  so  that, 
for  any  ray  p,  $(^,p)  can  be  computed  in  0{\ogn)  time.  The  preprocessing  time  of  their 
algorithm  has  been  reduced  to  O(nloglogn)  by  Guibas  et  al.  [GHLST].  For  the  general 
case,  however,  the  ray  shooting  and  other  visibility  problems  are  much  harder  even  for 
non-intersecting  segments.  For  example,  a  result  of  Suri  and  O'Rourke  [SO]  shows  that 
the  portion  of  a  non-simple  polygon  visible  from  a  fixed  edge  can  have  f2(n'')  edges  on  its 
boundary,  while  for  simple  polygons  such  a  region  is  bounded  by  only  0{n)  edges. 

We  axe  not  aware  of  any  ray-shooting  algorithm  for  non-simple  polygons  (or  for  an 
arrangement  of  segments),  which  answers  a  query  in  (9(log°'^'  n)  time,  using  roughly  linear 
space.  If  we  allow  quadratic  space,  then  a  query  is  easy  to  answer  in  time  O(logn)  (see 
Section  4.1).  Our  goal  in  this  paper  is  to  obtain  efficient  solutions  that  use  roughly  linear 
space,  and  to  establish  a  trade-off  between  space  and  query  time. 

For  a  special  case,  where  ^  is  a  set  of  fines,  a  result  of  Edelsbrunner  et  al.  [EGH*]  impfies 
that  one  can  construct,  in  randomized  expected  time  0{n'^^^  log'^  n),  a  data  structure  of  size 
O(nlog^n),  so  that  a  ray  shooting  query  in  A{C)  can  be  answered  in  0(>/nlog^n)  time. 
(The  preprocessing  has  been  made  deterministic  and  the  query  time  has  been  reduced 
to  0{y/n log n)  in  [Age].)  Unfortunately,  this  algorithm  does  not  apply  to  segments.  An 
algorithm  with  a  subfinear  query  time  for  the  case  of  segments  can  be  developed  using  the 
"recursive  space-cutting  tree"  of  Dobkin  and  Edelsbrunner  [DE]  (see  also  [EW]).  The  best 
known  algorithm  for  computing  $(^,  p)  is  by  Guibas  et  al.  [GOS],  which  constructs  a  data 
structure  of  size  0{n),  so  that  a  query  can  be  answered  in  0(n^/^+*)  time,  for  any  6  >  0. 
Their  algorithm  is  based  on  the  random  sampling  technique  of  [01],  [HW],  and  constructs 
a  multi-level  partition  tree.  The  preprocessing  of  their  algorithm  is  randomized  with 
0(n  log  n)  expected  running  time.  However,  the  preprocessing  can  be  made  deterministic 
without  any  additional  overhead  using  the  recent  algorithms  of  Matousek  [Maa]  or  Agarwal 
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[Agb]. 

In  this  paper  we  show  that  ray  shooting  can  be  performed  in  roughly  (that  is,  up  to 
poly-log  factors)  y/n  time,  while  still  using  only  roughly  linear  space  and  employing  deter- 
ministic, rather  than  randomized,  preprocessing  techniques.  We  first  give  an  algorithm  for 
the  case  of  non-intersecting  segments.  This  algorithm  constructs,  in  time  0(n^/^  log'^  n), 
a  data  structure  of  size  O(nlog^n)  so  that,  for  a  given  ray  p,  ^{Q,p)  can  be  computed  in 
0{y/n\og^  n)  time,  where  w  is  a  constant  <  4.33.  Our  algorithm  is  simpler  than  that  of 
[GOSl  because  it  maintains  only  a  two-level  data  structure.  We  then  extend  the  above  algo- 
rithm to  general  arrangements  of  segments.  Although  the  basic  idea  remains  the  same,  we 
need  several  new  techniques,  and  the  algorithm  is  more  complex.  In  this  case  a  query  can 
be  answered  in  0{yjna{n)  log^  n)  time,  using  0{na{n)\og*  n)  space,  after  ©(n^-'Mog'^  n) 
preprocessing.  Another  major  difference  between  the  two  cases  is  that  in  the  first  case 
we  can  report  all  K  intersections  between  a  query  ray  p  and  Q  in  0{y/n\og^  n  +  Klogn) 
time  while  we  still  do  not  know  how  to  report  these  intersections  in  a  comparably  efficient 
manner  in  the  general  case.  One  disadvantage  of  our  edgorithms  over  those  of  [GOS],  [DE] 
is  that  our  preprocessing  time  is  roughly  n^/^  instead  of  roughly  linear.  This  is  the  price 
that  we  have  to  pay  to  achieve  deterministic  preprocessing  and  to  reduce  the  query  time. 

The  second  problem  for  which  we  give  an  efficient  algorithm  using  the  spanning  tree 
data  structure  is  implicit  point  location.  The  implicit  point  location  problem  is  an  extension 
of  the  widely  studied  plzmax  point  location  problem  (see  [Ki],  [ECS],  [ST]).  In  the  latter 
problem,  a  planax  map  M  consisting  of  n  faces  is  given,  and  the  goal  is  to  preprocess  M 
into  a  data  structure  that  supports  fast  point  location  queries,  i.e.  queries  that  seek  the  face 
of  M  containing  a  query  point  p  .  The  above  algorithms  construct,  in  time  0(n  log  n)  (or 
sometimes  linear),  a  data  structure  of  line2Lr  size,  so  that  a  query  point  can  be  located  in  M 
in  O(logn)  time.  In  the  implicit  point  location  problem  the  map  is  defined  implicitly.  In 
particular,  we  assume  that  it  is  defined  as  the  axramgement  (i.e.  overlay)  of  a  given  set  of  n 
geometric  polygonal  (possibly  intersecting)  objects  of  some  simple  shape  (or  as  a  collection 
of  arbitrary  line  segments),  and  the  goal  is  to  obtain  certain  information  related  to  the 
arrangement  of  the  objects;  for  example,  to  determine  whether  a  query  point  lies  in  the 
union  of  the  objects.  A  more  formal  description  is  given  in  Section  7.  Guibas  et  al.  [GOS] 
have  presented  an  algorithm  with  0(n^/^'*"*)  query  time,  for  any  6  >  0,  using  the  rsmdom 
samphng  technique.  We  improve  the  query  time  to  0(\/n  \og^  n)  and  use  deterministic 
preprocessing.  The  algorithm  of  [GOS]  uses  0{n)  space,  while  ours  requires  O(nlog^n) 
space. 

Guibas  et  al.  [GOS]  have  described  several  applications  of  the  implicit  point  location 
problem,  such  as  polygon  containment,  implicit  hidden  surface  removal,  polygon  place- 
ment, etc.  We  show  that  our  implicit  point  location  algorithm  improves  the  query  time  of 
these  algorithms  as  well. 

This  paper  is  organized  as  follows.  In  Section  2  we  discuss  spemning  trees  with  low 
stabbing  numbers.  Section  3  describes  our  ray  shooting  algorithm  for  arrangements  of  non- 
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intersecting  segments.  In  Section  4  we  show  that  ray  shooting  queries  can  be  perforrned 
faster,  if  we  are  allowed  to  use  more  space.  Section  5  extends  the  algorithms  of  Sections 
3  and  4  to  report  all  intersections  between  Q  and  a  query  ray  p  at  logarithmic  cost  per 
intersection.  In  Section  6  we  generalize  our  ray  shooting  algorithms  to  arrangements  of 
arbitrary  (possibly  intersecting)  segments.  Section  7  presents  an  efficient  algorithm  for 
implicit  point  location  and  Section  8  discusses  other  applications  of  the  spanning  tree  data 
structure.  We  conclude  in  Section  9  with  some  final  remarks. 

2      Spanning  Trees  of  Low  Stabbing  Number 

Let  S  be  a  set  of  n  points  in  R'',  and  let  T  be  a  spanning  tree  of  S  whose  edges  are  line 
segments.  The  stabbing  number  (t{T)  of  T  is  the  ma.ximum  number  of  edges  of  T  that  can 
be  intersected  by  a  hyperplane  h.  Chazelle  and  Welzl  [CW]  (see  also  [We])  have  proved 
that,  for  any  set  of  n  points  in  IR  ,  there  exists  a  spanning  tree  with  stabbing  number 
0(n^"^^''),  and  that  this  bound  is  tight  in  the  worst  case.  For  a  family  T  of  trees,  the 
stabbing  number  cr{T)  is  5  if  for  each  hyperplane  h  there  is  a  tree  T  £  T  such  that  h 
intersects  at  most  s  edges  of  T.  . 

Chazelle  and  Welzl  [CW]  also  proved  that  a  spanning  tree  of  n  points  in  R''  with 
stabbing  number  0{n^~^^'')  can  be  constructed  in  polynomial  time.  In  the  plane,  a  span- 
ning tree  with  stabbing  number  0{^/n)  can  be  constructed  in  0(n'' log n)  time.  A  recent 
algorithm  of  Matousek  [Mab]  improves  the  running  time  to  0(n^/'^  log^  n)  at  the  cost  of 
increasing  a(T)  to  0{^/nlogn).  As  for  constructing  a  family  of  spainning  trees,  Edels- 
brunner  et  al.  [EGH*]  have  presented  a  randomized  algorithm,  with  expected  running 
time  O(n^/^log^n),  to  compute  a  family  T  =  {Ti,  . . .  ,Tk}  of  O(logn)  spanning  trees, 
with  cr(T)  =  0{y/n  \og^  n).  The  running  time  of  their  algorithm  has  been  improved  to 
0(n'*/^log^n)  in  another  rcindomized  algorithm  by  Matousek  [Mab].  (The  stabbing  num- 
ber of  T  computed  by  Matousek 's  algorithm  can  actually  be  improved  to  0{y/n\ogn)-  see 
[Age].)  An  additional  property  of  the  algorithms  of  [EGH*]  and  [Mab]  is  that  the  trees 
they  produce  are  actually  spanning  paths. 

The  best  known  deterministic  algorithm  for  constructing  a  family  of  spanning  path  is 
due  to  Agarwal  [Age],  who  has  shown  that 

Theorem  2.1  ([Age])  Given  a  set  S  of  n  points  m  the  plane,  one  can  deterministically 
construct  a  family  C  o/0(logn)  spanning  paths  on  S  witha{C)  =  0{s/n),  in  O(n^/Mog'*' n) 
time,  using  0{n^^'^)  working  storage,  where  u  is  a  constant  <  4.33.  Moreover,  for  any  query 
line  i,  one  can  determine  m  O(logn)  time  a  spanning  path  C  G  C  such  that  i  intersects  at 
most  0{^yn)  edges  of  C. 

D 

The  paths  constructed  by  [Mab]  and  [Age]  can  generally  be  self-intersecting.  However, 
Edelsbrunner  et  al.  [EGH*]  have  shown  that  a  spanning  tree  T  can  be  converted  into  a  non 
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self-intersecting  polygonal  path  C  in  O(nlogn)  time,  so  that  if  a  line  £  intersects  T  at  5 
points,  then  it  intersects  C  in  at  most  2^  points.  Therefore,  if  desired,  we  can  assume  that 
the  spanning  paths  produced  by  the  techniques  of  [Mab],  [Age]  are  non  self-intersecting. 


Figure  2:  C  and  B{C):  black  nodes  of  B  denote  Vb{C)  » 


Let  C  be  a  spcinning  path  on  5.  For  our  applications  we  need  to  convert  C  into  a 
balanced  binary  tree  whose  leaves  store  the  points  of  S  in  their  order  along  C.  Each  node 
V  oi  B  IS  associated  with  the  subpath  C„  of  C  connecting  the  points  stored  at  the  leaves  of 
the  subtree  rooted  at  v;  let  us  denote  by  Sv  the  subset  of  5  consisting  of  these  points  (see 
figure  2). 

A  line  £  stabs  a  node  u  of  5  if  ^  intersects  C„.  Let  Vb{£)  denote  the  set  of  nodes  v  oi  B 
such  that  V  is  not  stabbed  by  £  but  its  parent  (if  one  exists)  is  stabbed.  It  is  easily  seen 
that  {5v  :  V  6  Vb{£)}  is  a  disjoint  partitioning  of  S.  Moreover,  one  easily  verifies 

Lemma  2.2  //  a  line  £  intersects  s  edges  of  C,  then  \Vb{£)\  <  25logn  and  the  nodes  of 
Vb{£)  he  on  at  most  2s  paths  of  B. 


Another  simple  but  key  observation  is  that 

Lemma  2.3  A  line  £  intersects  a  polygonal  path  C  if  and  only  if  £  intersects  the  convex 
hull  of  the  vertices  of  C. 
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Lemma  2.3  implies  that  i  stabs  a  node  v  if  and  only  if  C  intersects  the  convex  hull  of  Sy. 
Since  an  intersection  between  a  Hne  and  a  convex  n-gon  can  be  detected  in  O(logn)  time, 
it  follows  that  Vb{^)  and  the  paths  containing  its  nodes  can  be  computed  in  0{  |V0(^)|  log  n) 
time,  if  we  store  the  convex  hull  of  the  subpath  Cy  at  each  node  v  of  B.  In  the  next  section 
we  will  show  that  this  computation  can  actually  be  done  in  time  O(|V0(^)|  +  logn),  using 
fractional  cascading  (cf.  [CGc]). 

3      Ray  Shooting  in  Arrangements  of  Non-intersecting 
Segments 

In  this  section  we  present  an  algorithm  that  preprocesses  a  given  set  ^  of  n  non-intersecting 
segments  so  that,  given  a  query  ray  p  emanating  from  a  point  p  in  direction  d,  ^{Q,  p)  can 
be  computed  quickly.  (For  technical  reasons  we  consider  p  as  an  open  ray,  i.e.  the  point  p 
does  not  belong  to  p.)  We  will  also  use  $(^,p)  to  denote  the  distance  of  that  point  from 
p:  if  no  such  intersection  exists,  we  put  $(^,p)  =  -|-oo.  Without  loss  of  generality,  we 
restrict  our  attention  to  rightward-directed  rays;  leftward-directed  rays  can  be  handled  in 
a  symmetric  way.  Denote  the  set  of  endpoints  of  the  segments  of  ^  as  S  =  {pi,  ...  ,Pm}i 
where  m  <  2n.  For  simplicity  of  exposition,  we  assume  that  no  segment  is  vertical.  Let  C 
=  {Ci ,  ...  Xk}  denote  a  family  oi  k  =  O(log  n)  spanning  paths  on  5,  with  a( C)  =  0(  v^)- 
We  show  how  to  preprocess  a  single  path  C  G  C.  First,  construct  B  =  B{C).  A  segment 
e  is  associated  with  a  node  u  of  B  if  one  of  its  endpoints  is  in  Sy.  Let  Ql  (resp.  QD  denote 
the  set  of  those  segments  of  Q  whose  left  (resp.  right)  endpoint  is  in  Sy  (see  figure  3;  note 
that  these  two  sets  are  not  necessarily  disjoint).  Let  i  denote  the  line  containing  the  query 
ray  p.  Since  |J  {QluQl)  =  Q,-we  have 
veVB(o 

$(a,p)   =     mm{mm{^{Ql,p),^gi,p)}}.  (3.1) 

Note  that  since  Cy  is  a  connected  path,  either  all  points  in  Sy,  for  a  node  v  e  Vb(^), 
lie  above  i  or  all  of  them  he  below  £,  say  above.  We  will  show  below  that  ^{G'y,  p)  and 
^(^v'P)i  for  V  e  Vb(£),  can  be  computed  in  O(logn)  time.  First,  a  few  notations:  Let 
i~  (resp.  i^)  denote  the  half  plane  lying  below  (resp.  above)  the  hne  L  We  distinguish 
between  the  two  sides  of  a  segment  e,  the  top  (resp.  bottom)  side  of  e  is  denoted  by  e+ 
(resp.  e").  We  say  that  a  ray  p  hits  e  from  above  (resp.  below)  if  slightly  to  the  left  of  their 
intersection,  p  lies  above  (resp.  below)  e.  If  we  think  of  e  as  expanded  into  a  very  thin 
rectangle  and  of  e+,  e"  as  denoting  the  top  and  bottom  sides  of  that  rectangle,  respectively, 
then  p  hits  e  from  above  if,  when  traversed  from  left  to  right,  p  first  intersects  e"*"  and  then 
e~,  and  symmetrically  for  rays  that  hit  e  from  below  (see  figure  4).  If  p  hits  e  from  above 
(resp.  below),  then  we  also  say  that  it  hits  e+  (resp.  e~). 

Ray  shooting  &  other  applications  May  11,  1989 


Ray  Shooting  in  Arrangements  of  Non-intersecting  Segments 


Figure  3:  Ql  and  Q^:  dcished  path  denotes  C^;  bold  lines  denote  Ql;  bullets  denote  Sy 


Figure  4:  Two  sided  segments:  p  hits  Cj  from  above  and  62  from  below 
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Since  the  segments  of  Ql,  for  v  €  Ve{£),  have  the  property  that  their  left  endpoints  He 
in  the  half  plane  f"*",  it  can  be  shown  that 

Lemma  3.1   Let  v  be  a  node  ofVei^)-   For  any  segment  e  G  G[,  P  either  does  not  hit  e.  or 
hits  e  from  below. 

D 

Before  proceeding,  we  introduce  a  linear  ordering  among  the  segments  of  Q^,  as  defined 
in  [GOS]  (see  also  [GY]).  As  we  will  see  later,  this  ordering  sorts  the  segments  in  a  manner 
that  is  consistent  with  any  order  in  which  they  can  be  crossed  by  a  rightward-directed  ray 
(from  below). 

Definition   3.2  [GOS]:     For  a  given  set  ^  =  {ci,  .  . .  ,e„}  of  segments, 

(i)  e,  <  Cj  if  there  exists  a  (non- vertical)  line  i  hitting  both  e~  and  ej  such  that  its 

s 

intersection  with  e,  lies  to  the  left  of  its  intersection  with  Cj ,  and  such  that  ^  does 
not  hit  any  ejj",  for  k  ^  i,j,  a,t  a.  point  between  e,  and  Cj. 

(ii)   e,<tj  if  there  exists  a  vertical  line  intersecting  both  e,  and  Cj  such  that  its  intersec- 
tion  with  e,  lies  below  its  intersection  with  Cj. 

(iii)  e,  <  Cj  if  e,  and  €_,  have  non-overlapping  x-projections  and  the  projection  of  e,  lies  to 
the  left  of  that  of  e^. 

(iv)  e,  <  Cj  (aJso  denoted  simply  as  e,   <  e^)  if  either  e,  precedes  e^  in  the  transitive 
closure  <  of  <,  or  e,  and  Cj  are  not  related  by  <  and  e,  <  e,. 

Theorem  3.3  ([GOS])    <  is  a  partial  order,   and    <    is  a  linear  order  that  extends  <. 
Moreover    <    can  be  computed  in  time  O(nlogn). 

n 

Remark  3.4:  It  is  possible,  for  a  pair  of  segments  Cj,  62  that  62  <  ei  within  some  set  Q, 

v*/z 

but  ei   <  62  relative  to  a  subset  Q'  C  G  (see  figure  5).  Therefore  it  is  important  to  mention 
the  set  relative  to  which  we  are  ordering  the  segments. 

Next  we  prove  a  technical  lemma  about    <    that  we  will  need  later.    Let  /^  (resp.  r^) 

v*/x 

denote  the  left  (resp.  right)  endpoint  of  a  nonvertical  segment  e. 

Lemma  3.5   For  all  segments  a,b  e  Gl  ,  if  r^  lies  below  £  and  x(ra)  <  x{lb),  then  a  <  b 
(relative  to  GD- 
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ei 


ei    <   62 


62    <    ei 


Figure  5:  Ordering  of  a  pair  of  segments  is  relative  to  a  set 

Proof:  Suppose,  to  the  contrary  that,  there  is  a  pair  of  segments  a,b  ^  Ql  such  that  r^ 
Hes  below  £  and  xir^)  <  x(li,),  but  b  <  a.  If  a  and  b  are  not  related  by  <,  then  a  <  6.  bv 

definition,  therefore  b  <  a  implies  that  b<a.  Thus  there  exists  a  sequence  of  segments  in 
Q\^  such  that  6  =  ei  <  62  <  •  •  •  <  e^  =  a.  Let  TT^^a  denote  a  shortest  sequence  axnong  all  such 

V  V  V 

sequences,  and  let  d^^a  denote  the  length  of  iTb.a-    We  obtain  a  contradiction  by  showing 
that,  for  every  A:  >  0,  there  is  no  sequence  Trta  such  that  (f(,,a  =  k. 


Figure  6:  Illustration  for  Lemma  3.5  db,a  =  3 

Obviously  db,a  >  2,  because  x(ra)  <  x(lb).  If  db,a  =  3,  then  there  is  a  segment  c  E  Ql  such 
that  b<c<a.  This  impHes  that  x{rc)  >  x{lb)  >  x{ra)  >  x(/c).  Let  q  be  the  intersection 
point  of  c  and  the  vertical  Une  x  =  x{ra).  Note  that  a  and  c  satisfy  the  following  properties: 
(i)  c<a,  (ii)  ar(rc)  >  x(ra),  (iii)  a  does  not  intersect  c,  and  (iv)  the  point  Ic  Hes  above  i 

(because  c  G  Si).  Using  these  properties  it  can  be  easily  proved  that  the  segment  qr^  hes 
below  £  (see  figure  6),  which  contradicts  the  assumption  that  b<c  (because  x{rc)  >  x{lb) 

V 

and  c  lies  below  b  at  x  =  x{lb)).  Hence  db,a  >  3. 

Now  assume  that,  for  all  segments  a,b  e  Q[  satisfying  the  properties  of  the  lemma, 
either  a  <  6  or  db^a  >  k.  Suppose  there  exists  a  pair  a,  b  such  that  b<a  and  db,a  =  ^-  Let 

b  =  ei  <  •  •  •  <  Cfc-i  <  Cyt  =  a  be  a  corresponding  shortest  sequence  7r(,,a,  and  let  c  =  Ck-i- 
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Since  7r(,,a  is  a  shortest  sequence,  it  is  easily  seen  that  x{rc)  >  x{ra).  Indeed,  let  Cj  be 
the  first  segment  in  this  sequence  whose  x-projection  overlaps  that  of  a.  Then  e^  must 
lie  below  a,  for  otherwise  we  would  have  obtained  a  cycle  in  <.  Hence  Cj  <  a  and  we  can 
shortcut  the  sequence  after  Cj.  Clearly,  x(iVj)  >  x{ra)-  Let  q  be  the  intersection  of  c  with 
X  =  x(ra),  as  above.  Again  we  can  argue  that  qT\  lies  below  (.  If  x(rc)  <  x(/i,),  then  c  and 
b  satisfy  the  property  of  the  lemma,  and  thus  contradict  the  inductive  hypothesis  because 
db  c  <  k.  On  the  other  hand  if  x(r<:)  >  x(/6),  then  we  have  c<b  (because  c  lies  below  b  at 
X  =  x(/(,)),  contradicting  the  assumption  that  b<c. 

Hence,  we  can  conclude  that  a  <  b. 

D 

Using  Lemma  3.1  and  Theorem  3.3  we  obtain 

Lemma  3.6   Let  £  =  (ci,  ...  .em)  denote  the  segments  of  Q'   ordered  with  respect  to    < 

(relative  to  Ql),  and  suppose  <^{Ql,p)  =  pPtej,  for  some  I  <  f  <  rn.  Then  for  all  i  <  f, 
e,  does  not  intersect  p. 

Proof:  If  e,  intersects  p,  it  does  so  at  a  point  to  the  right  of  p  fl  e/.  This  implies  that 
Cf  <  e,,  because  p  cannot  hit  any  e^  €  Ql  from  above,  so,  by  Theorem  3.3,  e/   <  e,,  or 

/  <  2,  a  contradiction. 

D 

Hence  the  original  problem  is  reduced  to  the  following  restricted  problem: 

Given  a  sequence  £  of  m  segments  sorted  according  to    <  ,  preprocess  £  so 

that  for  any  (rightward-directed)  query  ray  p  emanating  from  a  point  p  and 
lying  on  a  line  that  passes  below  the  left  endpoints  of  all  segments  in  £,  we  can 
quickly  determine  the  hrst  segment  ej{p)  of  £  hit  by  the  ray  p. 

A  possible  approach  to  solving  this  problem  is  to  do  a  binary  search  on  5,  where  each 
step  of  the  secirch  tests  whether  p  intersects  a  segment  in  some  contiguous  subsequence  of 
t  segments  of  £.  If  p  were  a  full  line  £,  then  such  an  intersection  could  be  easily  detected 
in  O(logi)  time  after  0{t\ogt)  preprocessing  (in  which  we  construct  the  convex  hull  of 
the  right  endpoints  of  these  t  segments).  However,  no  equally  fast  procedure  is  known  to 
detect  an  intersection  between  a  ray  and  such  a  set  of  segments.  To  overcome  this  problem, 
we  next  show  how  to  reduce  the  intersection  detection  problem  to  one  involving  the  hne 
containing  p  rather  than  p  itself. 

For  any  point  q  in  the  plcme,  let  e^  =  e/,(,)  denote  the  first  segment  of  £  whose  left 
endpoint  hes  to  the  right  of  (or  above)  q  (see  figure  7),  and  let  e^  =  £„(,)  denote  the 
segment  in  £  lying  immediately  above  g,  that  is,  the  vertical  ray  emanating  from  q  in  the 

Ray  shooting  &  other  applications  May  17,  1989 


Ray  Shooting  in  Arrangements  of  Non-intersecting  Segments  11 


Figure  7:  Segments  e/i,  e^  and  e^^ 

upward  direction  hits  e^  before  ciny  other  segment.  If  e^  (resp.  e^)  is  not  defined,  we  put 
h  =  m  +  \  (resp.  u  =  m  -\-  \).  Finally,  put  (^,  =  min{/i,u}. 

To  compute  e/j,  construct  a  balanced  binaxy  tree  L  whose  leaves  store  the  segments  of 
£  in  their  order  in  £.  For  each  interior  node  2  of  £-  we  store  the  rightmost  left  endpoint  of 
the  segments  stored  at  the  subtree  rooted  aX  z.  L  can  be  constructed  in  0(m  log  m)  time, 
and  it  is  easily  seen  that  th  caji  be  determined,  by  searching  for  q  through  L,  in  O(logm) 
time.  As  for  e„,  we  can  easily  calculate  it  in  time  O(log  m)  after  0(m  log  m)  preprocessing, 
as  in  [ST]. 

Lemma  3.7  The  query  ray  p  emanating  from  a  'point  p  cannot  intersect  any  segment 
e,  G  £,  for  i  <  4>p.  Moreover,  p  intersects  e,  for  i  >  (pp  if  and  only  if  its  right  endpoint  lies 
below  the  line  containing  p. 

Proof:  If  the  first  part  of  the  lemma  were  not  true,  then  there  would  exist  a  segment  e,, 
for  i  <  <t>p,  intersecting  the  ray  p.  In  this  case  the  left  endpoint  of  e,  must  lie  to  the  left 
of  p,  so  the  vertical  ray  77  from  p  in  the  upward  direction  must  intersect  e,.  But  then  the 
first  segment  e^  hit  by  rj  must  satisfy  i  >  u  >  (f>p  (because  e^  <  e,  and  by  definition  oi  4>p), 
a  contradiction  that  proves  the  first  half  of  the  lemma. 

The  "only  if"  peirt  of  the  second  half  of  the  lemma  follows  from  the  fact  that  if  both 
the  left  ajid  the  right  endpoints  of  a  segment  e  he  above  £,  then  e  cannot  intersect  £.  For 
the  "if"  part  let  e,  be  a  segment  of  £,  for  i  >  4>p,  whose  right  endpoint  hes  below  £,  but 
e,  does  not  intersect  p.  If  the  left  endpoint  /e.  of  e,  lies  to  the  right  of  p,  then  obviously  e, 
intersects  p,  so  Z^,  must  lie  to  the  left  of  p.  Since  e,  does  not  intersect  p,  the  intersection 
point  ^  of  e,-  and  i  lies  to  the  left  of  p.  Moreover  if  2;(re, )  <  x{U^^^^),  then  by  Lemma  3.5 
^i  <  Su(p)-  If^  3:(reJ  >  i(/e  ),  then  e,  and  e^ip)  must  have  x-projections  that  overlap  at 
some  point  between  ^  and  p\  since  Cj  lies  below  eu(p)  at  this  point,  we  again  have  e,  <  e^(p). 

V* 

Similarly  we  can  show  that  e,  <  e/,(p).  Hence  i  <  min{u,  h],  contradicting  the  eissumption 
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Figure  8:  Convex  hull  CE(Sp)  intersecting  £' 


that  i  >  4> 


V 


Lemma  3.7  implies  that  the  binaxy  search  technique  proposed  above  will  work,  provided 
we  can  detect  quickly  whether  any  right  endpoint  of  a  segment  in  some  subsequence  of  £ 
lies  below  L  In  other  words,  the  problem  now  has  been  reduced  to  that  of  detecting  an 
intersection  between  a  set  of  points  and  a  query  half  plane.  Clearly,  this  is  equivalent  to 
detecting  an  intersection  between  the  convex  hull  of  these  points  and  the  half  plane  (see 
figure  8). 

We  axe  now  ready  to  describe  how  to  preprocess  £  so  that  e/(p)  can  be  computed 
quickly,  for  any  ray  p  with  the  above  properties.  Let  r,  denote  the  right  endpoint  of  e,  G  £, 
and  let  i?  =  {ri,  . . .  ,rT„}.  We  construct  a  binary  tree  T  on  i?  in  the  Scime  way  as  we 
constructed  5,  i.e.  the  points  r,  axe  stored  in  the  leaves  of  T  in  order,  and  each  node  w  of 
T  is  associated  with  the  subsequence  Fiuj  of  R  containing  all  points  stored  at  the  leaves  of 
the  subtree  rooted  at  w. 

At  every  node  u>  of  T,  we  store  the  convex  hull  of  R^.  Using  T  we  can  easily  determine 
e/(p)  in  time  0{\og^  m):  We  first  find  4)^,  as  described  above,  in  O(logm)  time.  Then 
we  treat  the  suffix  {r^^,  ...  .r^]  of  i?  as  the  union  of  logm  subsets  R^,  w  G  T,  which 
we  can  compute  in  O(logm)  time.  We  test  each  R^  in  increasing,  left-to-right  order,  to 
find  the  first  w  for  which  the  line  £  containing  p  intersects  the  hull  of  il^,.  Then  we  do  a 
binary  search  within  R^  until  we  find  e/(p).  All  this  takes  O(log^m)  time.  However,  we 
can  easily  reduce  the  time  to  O(logm),  using  fractional  cascading.  This  is  possible  since, 
as  in  [CGc],  detecting  intersection  between  £  and  a  convex  polygon  axnounts  to  searching 
for  the  slope  of  £  in  the  sequence  of  slopes  of  the  edges  of  the  polygon  (see  [CGb],  [CGc] 
for  more  deteiils).  Therefore  we  have 
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Lemma  3.8  Given  a  set  £  of  m  non-intersecting  line  segments  in  the  ■plane,  we  can 
preprocess  them,  in  time  0(m  log  m),  into  a  data  structure  of  size  0(m  log  m)  so  that, 
given  a  (rightward- directed)  query  ray  p  whose  containing  line  lies  below  the  left  endpoints 
of  all  the  segments  m  S,  we  can  compute  the  first  segment  of  £  hit  by  p  m  time  O(iogm). 

D 

Remark  3.9:  Using  the  same  fractional  cascading  technique,  we  can  compute  Vb(^),  for 
any  line  £,  in  0(|Vb(£)|  +logn)  time  instead  of  O(|V0(^)|  logn)  time. 

Returning  to  the  original  problem,  Lemma  3.8  and  the  preceding  discussion  imply  that 
we  can  compute  $(^',,p)  (and  symmetrically  ^{Ql,p)),  for  each  u  G  Vs(£),  in  time  O(logn). 
Equation  (3.1)  and  Lemma  2.2  then  imply 

Theorem  3.10  Given  a  set  Q  of  n  non-intersecting  line  segments,  we  can  preprocess  it 
in  time  0(n^^^  \og'^  n),  for  some  u  <  4.33,  into  a  data  structure  of  size  0(n  log  n),  using 
0{n^^^)  working  storage,  so  that,  given  a  query  ray  p,  its  first  intersection  $(^,p)  with  Q 
can  be  computed  m  time  O(i/nlog   n). 

D 

Remark  3.11: 

(i)  The  space  used  can  be  reduced  to  0(n log'^n),  without  affecting  the  query  time,  if 
we  use  a  single  tree  structuire  instead  of  a  family  of  O(logn)  trees.  But  then  the 
preprocessing  time  increases  to  0(n"^ log n)  (see  [EGH*]). 

(ii)  If  we  allow  randomization,  the  (expected)  preprocessing  time  can  be  reduced  to 
0(n''/^log^n)  using  Matousek's  algorithm  [Mab],  but  then  the  query  time  bound 
increases  by  a  factor  of  logn. 

4     Tradeoff  between  Space  and  Query  Time 

In  this  section  we  show  that  the  query  time  for  the  ray  shooting  problem  in  arrangements 
of  non-intersecting  segments  can  be  improved  if  we  allow  ourselves  more  storage.  Similar 
tradeoffs  have  been  obtained  for  severed  related  problems,  such  as  computing  a  face  in  an 
arrangement  of  hues  [EGH*]  and  simplex  range  searching  [Age],  [Chd].  The  maun  result 
of  this  section  is  an  algorithm  for  computing  ^{Q,p)  with  0(--^  log^''^  -^  +  logn)  query 
time,  using  0{m)  space,  where  nlog'^n  <  m  <  n^. 
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Figure  9:  A  segment  e  and  its  duaJ  e' 

4.1      The  case  of  quadratic  storage 

First,  we  show  that  if  we  allow  0[n^)  space,  the  query  time  can  be  reduced  to  O(logn). 

Let  ^  =  {ci,  . . .  ,  Cn}  be  a  collection  of  n  non-intersecting  segments.  The  dual  of  a 
segment  e  =  a6  is  a  double  wedge  e'  formed  between  the  dual  lines  a',  h'  of  a,  b  respectively, 
and  not  containing  any  vertical  line  (see  figure  9).  Dualize  all  segments  e  6  ^,  obtaining  a 
set  Q'  of  n  double  wedges.  Let  £*  denote  the  set  of  lines  bounding  the  double  wedges  of 
Q'  (i.e.  the  duals  of  the  endpoints  of  segments  in  Q).  Let  A{C')  denote  the  arrangement  of 
£',  and  let  Wj  be  the  set  of  double  wedges  of  Q'  containing  the  face  /  G  A[C').  Standard 
duality  arguments  yield: 

Lemma  4.1  Let  p  he  a  point  lying  in  the  interior  of  a  face  f  of  A{C').  Then  p'  intersects 
each  segment  e  G  Wj  transversally  at  an  interior  point,  and  is  disjoint  from  any  other 
segm,ent  of  Q . 

D 

Lemma  4.2  //  the  segments  of  Q  are  non-intersecting,  then  for  all  points  p  m  the  face  f 
of  A{C'},  the  line  p'  intersects  the  segments  of  wj  m  the  same  order. 

Proof:  Suppose  there  are  two  points  x  and  y  in  a  face  /  such  that  the  lines  x*  and  y' 
intersect  the  segments  of  wj  in  two  different  orders.  Since  the  segments  in  Q  are  non- 
intersecting,  rotating  x'  towards  y'  (in  the  direction  that  avoids  a  vertical  orientation)  we 
must  reach  a  line  p*  that  either  contains  a  segment  of  wj,  or  passes  through  an  endpoint  of 
a  segment  of  wj.  (Note  that  this  claim  does  not  hold  if  the  segments  can  intersect.)  The 
dual  p  of  p'  is  a  point  that  Hes  on  the  segment  ly,  hence  in  /.  This,  however,  contradicts 
Lemma  4.1,  thus  showing  that  the  duals  of  all  points  in  /  intersect  the  segments  of  wj  in 
the  same  order. 
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D 

Sort  the  segments  in  wj  in  the  order  provided  by  Lemma  4.2.  For  a  ray  p,  let  the  image 
of  p  be  the  dual  of  the  line  containing  p.  If  the  image  of  a  ray  p  lies  in  the  face  /  €  A{C'), 
then  ^{Q,p)  can  be  computed  in  O(logn)  time  by  a  binary  search  on  Wf.  Therefore,  it 
suffices  to  show  how  to  store  all  the  lists  wj  using  only  0{n^)  space,  so  that  binary  search 
in  each  of  them  can  still  be  done  in  O(logn)  time. 

Let  V  denote  the  dual  graph  of  A{C'),  i.e.  the  graph  whose  nodes  represent  faces  of 
A{C')  and  whose  edges  connect  pairs  of  nodes  representing  adjacent  faces.  Let  T  denote 
a  spanning  tree  of  T>.  We  can  convert  T  into  a  path  11  by  tracing  an  Eulerian  tour  around 
the  tree.  Observe  that  if  two  vertices  ui,U2  in  11  represent  faces  /i,/2  sharing  an  edge  7, 
which  is  a  portion  of  a  line  £,  then  if/,  ©  w/^  is  the  set  of  segments  having  the  dual  of  ^  as 
an  endpoint.  Let  6^  denote  this  set  of  segments.  The  set  Wf^  can  be  obtained  from  wj^  by 
deleting  the  segments  of  6-,  D  wj^  and  inserting  the  segments  of  6^  —  Wf^ .  Therefore  we  can 
maintain  all  lists  wj  using  a  persistent  data  structure  (see  [Co],  [ST],  [DSST]).  Since  at 
each  edge  7  of  11,  only  the  segments  of  6^  axe  inserted  or  deleted,  the  totaJ  space  required  to 
store  all  wj  is  0(n  +  X^ygn  l^-cl)  ^"*i  the  totaJ  preprocessing  time  is  0((n  +  Xl7en  l<^-yl)log")- 
Moreover,  using  this  persistent  data  structure,  ^{Q,p)  can  be  computed  in  (9(logn)  time 
(see  [ST]).  Thus,  it  suffices  to  prove  that  J2-t  \^-y\  =  O(n^).  Suppose  the  segments  of  Q  have 
t  <  2n  distinct  endpoints  and  Ui  segments  are  incident  to  the  z"*  endpoint.  It  is  easy  to 
check  that  if  7  is  a  portion  of  the  line  dual  to  the  ?"'  endpoint,  then  \6-,\  <  Vi.  Obviously 
ELi  ^'  =  2n  and  each  hne  of  C  is  split  into  <t  +  1  edges,  which  implies  that 

El-^^l      <      D^  +  l)^.       <      2n(2n  +  l). 

-1&V  1=1 

Hence,  we  have 

Theorem  4.3  Given  a  collection  Q  of  n  non-interseciing  segments,  we  can  preprocess  Q, 
in  O(n^logn)  time,  into  a  data  structure  of  size  0{n'^)  so  that,  for  any  query  ray  p,  ^{Q,p) 
can  be  computed  in  O(logn)  tim.e. 

D 

4.2      The  general  case 

Theorem  3.10  and  Theorem  4.3  represent  roughly  two  extremes  of  the  spectrum,  because 
we  need  at  least  0{n)  space,  and  we  cannot  hope  to  answer  a  query  in  o(logn)  time.  The 
general  case  where  the  allowed  storage  m  assumes  an  intermediate  value  between  n  log^  n 
and  n^  is  handled  as  follows.  For  technical  reasons  we  assume  for  the  time  being  that  no 
endpoint  of  a  segment  in  Q  has  degree  >  3  (that  is  incident  to  more  than  three  segments  of 
Q).  In  Section  4.3  we  show  how  to  handle  degenerate  cases  (i.e.  when  there  are  endpoints 
of  degree  >  3). 
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Using  the  algorithm  of  [Agb],  partition  the  dual  plane,  in  time  (9(nr  log nlog'^"^  r),  into 
M  =  O(r^)  triangles  Ai,  .  .  .  ,  Am,  each  meeting  at  most  ^  lines  dual  to  the  endpoints  of 
the  given  segments,  where  r  is  a  parameter  to  be  chosen  later.    Let  £*  denote  the  set  of 

dual  lines  that  intersect  the  triangle  A,,  for  z  =  1 A/;  then  |£*|  <  7.    For  each  A,, 

define  the  subset  Q,  of  Q  to  consist  of  all  segments  e  having  at  least  one  endpoint  whose 
dual  is  in  £,•.  Obviously  |C^.|  <  7^.  We  also  define  W,  C  Q  as 

W,   =    {e|e  €  a  A  A,  C  e*}. 

It  is  easily  checked  that 

Lemma  4.4   For  each  point  p  lying  inside  A,,  the  line  p*  does  not  intersect  any  segment 

ofg-{g,uw,). 

D 

Lemma  4.4  implies  that 

$(^,p)   =   min{$(VF.,p),  $(^.,p)}  (4.1) 

where  A,  is  the  triangle  containing  the  image  of  p.  Using  the  same  argument  as  in 
Lemma  4.2,  we  can  prove 

Lemma  4.5  All  lines  whose  dual  points  lie  inside  A,  intersect  the  segments  of  W,  in  the 
sam.e  order. 

a 

We  can  thus  order  the  segments  of  Wi  in  the  order  provided  by  Lemma  4.5,  and  compute 
$(W,,  p),  for  amy  ray  p  whose  image  hes  in  A,,  in  O(logn)  time,  using  binary  search.  Let 
Ai  and  A2  be  two  adjacent  triangles  and  let 

Q12   =    {e|ee^i  A  A2  Ce*}.  (4.2) 

It  is  easily  seen  that  W2  =  {W^  U  Qn)  -  Qi-  Since  l^il,  1^21  =  O(^),  we  have  |H'2  -  l^'il 
=  O(^).  As  earher,  we  define  a  graph  D,  whose  vertices  are  the  triangles  A,  and  whose 
edges  connect  pairs  of  vertices  representing  adjacent  triangles.  Now  an  edge  between  Vi 
and  V2  has  the  set  Qx  U  Q2  associated  with  it.  Again,  we  construct  a  path  11  on  a  spanning 
tree  of  P,  and  obtain  a  persistent  data  structure  Ti(^)  to  store  W^  for  all  triangles.  It  can 
be  easily  shown  that  Ti(^)  requires  0{Tir)  space,  and  can  be  constructed  in  0(nr log n) 
time.  Now,  for  any  ray  p,  ^{W,,p)  can  still  be  computed  in  O(log  n)  time,  where  Ai  is  the 
triangle  containing  the  image  of  p. 

We  preprocess  each  ^,  into  a  data  structure  T2(^,)  of  size  0(|^,p)  for  ray  shooting 
queries,  as  described  in  Section  3,  so  that  for  any  ray  having  its  image  in  A,,  we  can  find 
$(^.,p)inO(yflog2^)time. 
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To  compute  ^{Q,p),  for  a  given  ray  p,  we  first  find  the  triangle  A,  that  contains  its 
image;  this  can  be  done  in  O(logn)  time,  using  an  efficient  point  location  algorithm  [EGS], 
[ST].  It  follows  from  (4.1)  that  ^{Q,p)  can  be  computed  by  calculating  each  of  $(Wi,p) 
and  ^{Q,,p),  as  described  above;  therefore  the  query  time  is 


Q{n)   =   0[J-\og'-+[ogn 


O  (,/-log^-  +  logn)  . 
VV  r  r  J 


As  for  the  space  complexity  S[tl)^  we  need  O^r^)  space  to  store  the  triangle  Ai,  ...  ,  A^, 
O^nr)  space  to  store  Ti  and  0(Mog"'  ^)  to  store  each  Q^  (cf.  Theorem  3.10).  Thus, 

5(n)   =    0(nr)  +  o('r2.-log3-)       =       ofnrlog^-V 


If  we  choose  r  = 


m 
— ; — ^ ,  then  S(n)  =  O^m)  and  the  query  time  becomes 

"log    7^ 


Q{n)      =      O 


=     O 


n  -      2     " 


"       1        7/2       "  1  \ 

log  '    —y=  +  log  n     . 


m 


Next,  we  bound  the  preprocessing  time  P{n).  We  can  compute  Ai ,  ...  ,  Am  in  0{nr  log  n- 
log'^"^  r)  time  (see  [Agb]).    Since  Ti  can  be  constructed  in  0(nr  log  n)  time  and  each  Q, 
can  be  preprocessed  in  O^^jf^log"^  ^)  time  (cf.  Theorem  3.10),  we  have 


P{n)      =      O(nrlognlog"'-^r)  +  0     r 


n 


3/2 


r-3/2 


log^ 


n 


=     O Tnr log n log"-' r  +  n^/2^ log"-] 


,_i  rn    ^     3^2 


+  n^ 


\   "log  7^ 
=     Ofm log"  n  +  nV^ log"-^/2"\ 


^  n  log    ^  v^ 


Since  we  need  0{nr)  space  to  compute  Ai,  ...  ,  Am,  and  O(^)  to  preprocess  each 
^,,  the  total  space  required  for  preprocessing  is 


O  \nr  + 


n 


3/2^ 


r-3/2 


=    o 


m 


n 


3/2 


l°«^fe       mVV(nV^)3/'^ 


=    o( 


n" 


V  "^  +       3/2 


log 


3/2      "     \ 

y/mj 


Hence,  we  can  conclude 
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Theorem  4.6  Given  a  collection  Q  of  n  non-intersecting  segments  in  the  plane  with  the 
property  that  no  endpomt  has  degree  >  3,  and  a  parameter  nlog^n  <  m  <  n'^,  we  can 
■pre-process  Q,  m  0(m\og'^  n  +  ns/mlo^'^^^  n)  time,  into  a  data  structure  of  size  0{m)  so 

that,  for  query  ray  p,  we  can  compute  ^(Q ,  p)  in  O     —r=\og'^'^—=  +  log  n      time.    The 

\  \/rn  \/rn  J 


working  storage  required  for  preprocessing  is  O  \^m  -\ —  log  '^  -y=  J   time. 


n"^    .      o/T     n 


Remark  4.7: 

(i)  If  we  allow  randomization,  then  using  Matousek's  algorithm  Q,  can  be  preprocessed  in 
O  (  (^]  log'^  n  J  time,  but  the  query  time  increases  by  a  factor  of  O(logn).  There- 
fore, following  the  same  ajialysis  we  obtain 

P(n)      =      O  (m^/^n^/^  +  m  log'^  n) 
Q{n)     =     of-^log^/^^  +  lognV 

(ii)   If  we  maintain  a  single  tree  data  structure  for  each  ^,,  the  query  time  can  be  reduced 
to  O(-^log^n),  but  the  preprocessing  time  increases  considerably. 

4.3      Coping  with  degenerate  cases 

The  analysis  of  the  algorithm  described  in  the  previous  subsection  breaics  down  if  the 
segments  of  Q  have  endpoints  of  axbitreirily  laxge  degree,  because  then  we  caxmot  guarantee 
that  \Q,\  =  O(^),  2Lnd  the  analysis  to  bound  the  totaJ  space  required  to  store  Ti  relies 
heavily  on  this  bound  for  \Q,\.  In  this  subsection  we  overcome  this  difficulty  by  showing 
that,  given  a  set  ^  of  n  non-intersecting  segments,  we  can  transform  it  into  another  set  Q' 
of  at  most  3n  (non-intersecting)  segments  such  that  no  endpoint  of  a  segment  in  Q'  has 
degree  >  3,  and  ^{Q,p)  can  be  determined  from  ^{Q',p)  in  0(1)  time. 

Let  ^p  =  {ci,  . . .  ,  Ct}  be  a  subset  of  segments  of  Q  all  having  a  common  endpoint  p.  Let 
6  be  the  minimum  distance  from  p  to  its  closest  neighbor  in  Q  —  Qp,  and  let  c  be  the  circle 
of  radius  |  with  p  as  center.  For  a  segment  e,  G  ^p,  let  q,  denote  the  intersection  point  of 
c  and  e,.  Assume  that  the  segments  of  Qp  axe  ordered  in  counter-clockwise  direction  along 
p.  There  are  two  cases  to  consider: 

(i)  There  exist  two  consecutive  segments  in  Qp,  say  e,  and  Cj,  such  that  the  angle  between 
et  and  Cj  is  >  180°.  For  1  <  z  <  f ,  we  remove  the  portion  of  e,  that  lies  in  the  interior 
of  c  (i.e.  pg7),  and  add  the  segments  q^,  . . .  ,qt-iqt  to  Q  (see  figure  10a). 
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(a) 


Figure  10:  Modifying  segments  having  a  common  endpoint  of  degree  >  3. 


(ii)  The  angle  between  every  two  consecutive  segments  of  Qp  is  <  180°.  For  each  e  G  Qp, 
we  remove  the  portion  of  e  that  Hes  in  the  interior  of  c,  and  add  the  segments 
qlgi",  . . .  ,qt-iqt,Mi  to  Q  (see  figure  10b). 

We  repeat  this  process  for  each  endpoint  of  the  segments  of  Q  whose  degree  is  >  3.  Let 
Q'  be  the  new  set  of  segments;  obviously  \Q'\  <  3n,  and  each  endpoint  has  degree  <  3.  It 
can  be  eeisily  proved  that,  for  a  ray  p,  if  ^{Q\p)  lies  on  a  segment  of  Q,  then  $(^,p)  = 
^{Q',p).  Otherwise  if  ^{Q',p)  Hes  on  a  segment  q,-iq,,  and  e,_i  (resp.  e,)  is  the  segment 
of  Q  incident  to  9,_i  (resp.  q,),  then  ^{Q,p)  lies  on  either  e,_i  or  e,.  Therefore,  ^{Q,p) 
can  be  computed  from  <E>(^',p)  in  0(1)  time.  Moreover  for  each  endpoint  p,  the  minimum 
distance  6p  can  be  computed  in  0(n  log  n)  time  by  constructing  the  closest  point  Voronoi 
diagram  of  Q  [Ya]. 

Hence  by  Theorem  4.6,  we  have 

Theorem  4.8   Given  a  set  Q  of  n  segments  and  a  parameter  nlog^n  <  m  <  n},  we  can 
preprocess  it,  in  time  0(777  log"^  n  +  n-y/mlog'^"^' ^ -^),  into  a  data  structure  of  size  0{m) 

so  that,  for  a  query  ray  p,  one  can  compute  ^{Q,p)  tn  0(-^log'^^^  -^  +  logn)  time. 

D 
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5      Reporting  All  Intersections 

In  the  last  two  sections  we  gave  algorithms  to  compute  ^{Q,p)  for  a  collection  of  non- 
intersecting  segments.  We  now  extend  these  algorithms  to  solve  the  following  problem: 

Given  a  set  Q  of  n  non-intersecting  segments,  preprocess  it  so  that,  for  a  query 
ray  p.  one  can  quickly  compute  all  intersections  Ip  between  Q  axid  p  in  their 
order  along  p. 

Dobkin  and  Edelsbrunner  [DE]  have  given  an  algorithm  that  preprocesses  Q  into  a  data 
structureof  hnear  size  so  that,  for  a  query  ray  p,  Ip  can  be  computed  in  0{n°-^^^  +  \Ip\)  time. 
(In  fact  their  algorithm  works  for  an  arbitrary  collection  of  segments.)  We  first  present 
an  algorithm  that  uses  roughly  linear  space,  by  generalizing  the  algorithm  described  in 
Section  3. 

Preprocess  Q,  as  in  Section  3,  in  O(n''^^log"' n)  time  using  O(nlog  n)  space.  For  a 
given  ray  p,  we  compute  Ip  as  follows.  Let  I  denote  the  line  containing  the  ray  p,  and  let 
C  be  a  spanning  path  in  C  that  intersects  £  in  0{\/n)  edges.  As  described  in  Section  2, 
compute  Vb(^)  in  0{\/nlogn)  time.  Now  we  report  all  intersection  points  in  Ip  by  wadking 
along  the  ray  p  and  stopping  at  each  point  of  Ip.  For  a  point  q  G  p,  let  p,  be  the  ray 
emanating  from  q  and  contained  in  p. 

The  algorithm  maintains  the  following  invariant:  When  we  are  at  a  point  q  E  p,  we 
maintain  a  hst  of  all  points  ^{Q^,,  Pg),  for  all  v  G  Vs(£),  as  a  priority  queue  Q  (with  respect 
to  their  order  along  p).  Observe  that  Q  remains  the  same  between  two  consecutive  points 
of  Ip,  and  that  the  root  of  Q  stores  the  point  of  Ip  that  we  are  going  to  encounter  next. 
Therefore,  it  suffices  to  show  how  to  update  Q  cifter  visiting  a  point  of  Ip.  Suppose,  when 
we  are  at  a  point  q,  the  root  of  Q  stores  a  =  <5(^u,p,),  for  some  u  6  Vsii).  It  is  easily 
seen  that  when  we  cross  a,  the  next  intersection  point  of  p  and  Q^,,  for  all  u  €  Vb(£)  -  {u}, 
does  not  change.  Thus  Q  can  be  updated  by  deleting  a  from  Q,  and  inserting  ^(Qu.Pa) 
in  Q  provided  <5(^u,p<,)  ^  oo.  Continue  this  process  until  Q  becomes  empty.  It  is  easily 
seen  that  this  procedure  reports  all  intersection  points  of  p  and  the  segments  of  Q  in  their 
order  along  p. 

To  bound  the  running  time  of  the  algorithm,  observe  that  initially  we  spend  0{y/n\ogn) 
time  to  construct  the  queue  Q  ior  q  =  p,  and  then  spend  O(logn)  time  in  updating  Q 
after  each  intersection.  Hence,  we  have 

Theorem  5.1  Given  a  collection  Q  of  n  non-mtersecting  segments,  we  can  preprocess  it, 
in  time  0(n^/^  log*^  n),  into  a  data  structure  of  size  O(nlog^n)  so  that,  given  a  query  ray 
p,  Ip  can  be  computed  m  0{y/n\og^n  +  \Ip\logn)  time. 

D 
An  immediate  corollary  of  Theorem  5.1  is 
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Corollary  5.2  Given  a  collection  Q  of  n  non-intersecting  segments,  we  can  preprocess 
it.  m  time  0{n'^'^  log^  n),  into  a  data  structure  of  size  0(nlog  n)  so  that,  given  a  query 
segment  e,  we  can  compute  all  K  intersections  between  e  and  Q  in  time  0{\/n\o^  n  + 
K  log  n ) . 

Next  we  show  that,  as  in  Section  4,  the  query  time  can  be  improved  if  we  allow  more 
space.  Now  preprocess  Q  as  described  in  Section  4  (if  the  segments  of  Q  have  endpoints 
with  degree  >  3,  we  modify  the  set  Q,  as  described  in  Section  4.3).  Recall  that,  in  Section  4, 
we  maintain  two  data  structures  —  (i)  the  persistent  data  structure  Ti  to  store  W,  for 
each  triangle  A,,  and  (ii)  ~T2{Gx)  for  ray  shooting  queries.  For  a  query  ray  p,  we  compute 
Xp  as  follows. 

Suppose  the  ray  origin  p  lies  in  the  triangle  A,,  let  the  sorted  W,  be  (ei,e2,  ...  ,6^) 
and  suppose  $(  W,,  p)  €  e^.  Then  by  Lemma  4.2,  e^,  ...  ,  e^  intersect  p  in  that  order  along 
p,  and  we  thus  obtain  all  intersections  between  \\\  and  p.  The  intersections  between  ^, 
and  p  cire  obtained  by  the  procedure  described  above,  except  that  the  the  size  of  Q  is  now 
only  0{J^\og  ^)  because  \Q,\  <  ^.  Ip  is  then  obtained  by  merging  the  two  output  Usts  of 
intersections  with  IV,  and  ^,.  Hence,  following  the  same  analysis  cis  in  Section  4,  we  can 
conclude 

Theorem  5.3   Given  a  collection  Q  of  non-intersecting  segments  and  a  parameter  nlog''  n  < 
m  <  n^,  we  can  preprocess  it,  in  time  0(mlog'^n  +  n^/rn\og'^~^/'^  -^),  into  a  data  struc- 
ture of  size  0{m)  so  that,  given  a  query  ray  p,  Jp  can  he  computed  in  0{-j=\o^^^  -^  + 
logn  +  |Jp|log^)  time. 

D 

Corollary  5.4    Given  a  collection  Q  ofn  non-intersecting  segments  and  a  parameter  n  log^  n 
<  m  <  n^ ,  we  can  preprocess  it,  in  time  C>(m  log*^  n  +  n-y/mlog"'""'''^  -^),  into  a  data  struc- 
ture of  size  0{m)  so  that,  given  a  segment  e,  we  can  compute  all  K  intersections  between 
€  andg  in  0{^\og'/'^  +  logn +  K  log  ^)  time. 

a 
6      Ray  Shooting  in  General  Arrangements  of  Segments 

In  this  section  we  extend  our  algorithm  to  genercd  arrangements  of  possibly  intersecting 
segments.  The  section  is  organized  as  follows.  In  Section  6.1  we  describe  how  to  preprocess 
Q  for  ray  shooting  queries,  and  in  Section  6.2  we  show  how  to  answer  a  query.  We  analyze 
the  time  and  space  complexity  of  our  algorithm  in  Section  6.3  and  finally  derive  a  tradeoff 
between  space  and  query  time,  similar  to  that  of  Section  4,  in  Section  6.4. 
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6.1      Preprocessing  the  segments 

In  this  section  Q  denotes  an  arbitrary  collection  of  n  segments  in  the  plane.  To  simplify 
the  exposition,  we  assume  that  the  segments  of  Q  are  bounded.  The  preprocessing  of  Q 
is  done  as  follows.  We  construct  a  partition  tree  T,  and  associate  with  each  node  v  E  T 
a  collection  Gv  Q  Q  oi  n^  segments,  a  triangle  A^,  and  another  auxiliary  set  G'y  of  n[ 
segments.  If  Uy  <  c.  for  some  fixed  constant  c,  then  u  is  a  leaf  of  T.  Otherwise  it  is  an 
internal  node  of  T,  which  is  further  processed  as  follows.  For  some  fixed  constant  r  >  2, 
partition  Ai,  into  M  =  0(r'^)  triangles  Ai,  ...  ,  Am,  using  the  algorithm  of  [Agb]  (or  of 
[Maa]),  so  that  each  triangle  A,  meets  at  most  ^  lines  containing  the  segments  of  Qy. 
Create  M  children  u'l,  ...  ,  u'^,/  of  u,  and  associate  with  each  child  w,  the  corresponding 
triangle  A^;,  =  A,.  We  put  a  segment  e  of  Qy  in  Q^,  if  at  least  one  of  the  endpoints  of  e 
lies  in  A,.  We  also  associate  with  tw,  an  auxiliary  set  Q'^  of  all  segments  of  Qy  that  cross 
A,  but  do  not  have  any  endpoint  inside  A,.  Let  Ali,  be  the  planar  map  formed  by  the 
triangles  Ai,  ...  .Am-  The  root  u  of  T  is  associated  with  Q  itself,  ajid  Au  is  a  triangle 
that  contains  all  the  segments  of  Q.  Moreover,  Q'^  =  0,  by  definition. 


(a):  Q.,  G'y  (b):  Ay,  By 

Figure  11:  (a):  solid  lines  denote  Q^;  (b):  solid  lines  denote  Ay 

We  preprocess  each  node  v  E  T  as  follows.  Preprocess  the  planar  map  My  for  point 
location  queries  (see  [EGS],  [ST])  and  store  the  resulting  data  structure  at  v.  Let  €[ 
denote  the  set  of  lines  containing  the  segments  of  ^(,;  \C'J  <  n'„.  Preprocess  £[  into  a 
data  structure  Ti{C'J  for  computing  ^{C'y,p),  using  the  algorithm  of  Edelsbrunner  et  al. 
[EGH*]  for  ray  shooting  in  arrangements  of  lines  (see  also  [Age]).  If  ^{C[,p)  Hes  outside 
Av,  then  we  reset  it  to  +oo. 


Next,  partition  the  segments  of  Qy  U  G'y  into  two  parts: 
(i)   Ay-,  "long"  segments  having  at  most  one  of  their  endpoints  inside  A^;  \Ay\  <  Uy  +  n 


(ii)  By-,    "short"  segments  having  both  of  their  endpoints  inside  A^,  that  is,  segments 
lying  fully  inside  A^;  since  |.4„|  >  n^,  we  have  \By\  <  n„. 

Definition  6.1:  The  zone  of  a  triangle  A  in  an  arrangement  AiQ)  of  a  set  Q  of  segments 
is  the  collection  of  the  face  portions  /  n  A,  for  all  faces  /  6  AiQ)  that  intersect  5 A  (see 
figure  12). 
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Figure  12:  Zone  of  a  triangle  A,;  some  faces  are  non-simple  polygons 

Using  the  same  argument  as  in  [EGP*],  it  can  be  proved  that  the  toted  number  of 
edges  in  the  zone  of  a  triangle  in  an  arrangement  of  n  segments  is  0{na{n)),  where  a{n) 
is  a  functional  inverse  of  Ackermajin's  function.  Let  H^  denote  the  zone  of  A^  in  A(A^). 
Since  \A^\  <  n^  +  n'^,  there  are  0{{ny  +  n'^)a{n^  +  n'„))  edges  in  Hy.  Observe  that  every 
segment  of  ,4^  intersects  dAy,  therefore  Hy  U5A„  is  a  collection  o{  k  <  ^  simple  polygons 
P  =  {7^1,  . . .  ,Vk},  each  of  which  touches  dAy  (see  figure  13).  In  other  words,  each  V, 
has  at  least  one  edge  that  is  a  portion  of  5Av.  Moreover,  a  point  p  £  dAy  —  UAy  lies  in 
exactly  one  polygon  "P  G  P.  Sort  the  intersection  points  between  segments  of  Ay  and  dAy 
so  that,  for  a  point  p  6  dAy  —  UAy,  the  polygon  containing  p  can  be  computed  quickly, 
using  binary  search.  Preprocess  each  simple  polygon  V  E  P  into  a  data  structure  for  ray 
shooting  queries  using  the  algorithm  of  Chazelle  and  Guibas  [CGa]  (see  also  [GHLST]).  If 
$CP,p)  lies  on  5At,  (that  is,  p  does  not  intersect  V  in  the  interior  of  A^),  then  ^{V,p)  is 
reset  to  +oo.  Let  T2(A„)  be  the  union  of  all  these  structures. 


Figure  13:  Zone  of  A„  for  A^,;  each  face  is  a  simple  polygon  (possibly  with  some  overlapping  edges) 

Finally,  we  preprocess  the  short  segments  in  By  for  ray  shooting  in  non-simple  polygons. 
Let  Ty  denote  the  boundary  of  the  unbounded  face  in  A{\Jj<m  By,  ),  where  Wi,  ...  ,  wm  are 
the  children  of  v.  By  the  result  of  Pollack  et  al.  [PSS]  (see  'also  [GSS]),  :F„  has  0{nya{ny)) 
segments,  as  \\J]<m  By,J  <  \By\  <  Uy.   Since  the  segments  of  j;  are  non-intersecting,  we 
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preprocess  them  into  a  data  structure  TsiBy)  for  computing  ^{J^y,p],  using  the  algorithm 
described  in  Section  3. 

We  repeat  this  preprocessing  for  every  node  v  of  T.  The  resulting  collection  of  data 
structures  is  the  output  of  the  preprocessing  stage. 

6.2      Answering  a  query 

Let  p  be  a  query  ray  emanating  from  a  point  p  in  direction  d.  The  query  is  answered  by 
traversing  a  path  lip  of  T  and  computing  a^  =  ^{Qv,  p)  at  each  node  t'  6  Hp  in  a  bottom-up 
fashion.  At  the  end  of  this  process  we  obtain  at  the  root  u,  ct^  =  ^{Qu,  p)  =  ^{G-p)- 

The  path  lip  is  defined  so  that  for  each  node  v  along  lip  the  ray  origin  p  lies  in  A^. 
The  processing  of  a  node  v  of  lip  is  done  as  follows.  If  f  is  a  leaf,  then  we  compute  a^ 
directly  using  a  brute-force  method;  otherwise  we  locate  the  triangle  A(  G  -M^  containing 
the  ray  origin  p  and  obtain  the  corresponding  child  t  of  v  along  lip.  We  next  test  whether 
^iQv,p)  lies  inside  A(.  This  is  done  by  calculating  each  of  cr^  =  ^{Qt,p}  and  a\  =  ^{Q[,p). 
The  former  of  the  two  is  computed  recursively  ajid  the  latter  is  computed  using  Ti{C[).  If 
$(£t,  p)  €  A(,  then  a[  =  ^{C[,  p),  because  segments  of  Q'^  do  not  have  their  endpoints  inside 
At;  moreover  ^{C[,p)  ^  A<  implies  a[  0  A(.  Thus  a[  can  be  determined  by  computing 
$(£;,p),  using  Ti(£;). 

If  either  at  or  cr[  lies  inside  A(,  then  we  have  already  found  a^  =  $(^y,p).  Otherwise 
we  know  that  a^  lies  outside  A(,  in  which  ca.se  we  proceed  a^  follows.  Let  Wp  denote  the 
set  of  children  of  v,  other  than  A(,  that  intersect  p.  Wp  can  be  computed  by  traversing 
the  ray  p  in  the  planar  map  Mv-  It  is  easily  seen  that 

$(a„p)      =      min{$(au,UCp)} 

=     min  {min{$(A^,/)),  $(5^,  p)]] 

=     °^n{mm{*(.4.,p)},  min  {$(5.,p)}} 

=     min|rnin{$(A^,p)},  {$(   IJ    B^^p)]\  (6.1) 

For  w  G  Wp,  p  does  not  lie  in  A^i,,  therefore  either  p  does  not  intersect  A^,  or  if  it 
does,  then  <I'(A^,p)  lies  on  jm  edge  of  H^.  Let  ^  be  the  first  intersection  point  of  p  and 
5Au,,  emd  let  V  be  the  polygon  of  H,^  containing  ^  (see  figiire  14).  (Assume  that  (,  ^  \JA^, 
because  then  we  already  know  $(A^,p).)  If  p'  is  the  ray  emanating  from  ^  in  the  same 
direction  as  p,  then  it  is  obvious  that 

$(A^,p)  =  ^{A^,p')  =  <5(P,p'). 

But  ^{'P,p'),  and  thus  also  $(Au,,p),  can  be  computed  from  T2{Au,)  in  the  logarithmic 
time. 
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Figure  14:  A{A^^)  and  V 


Since  p  does  not  intersect  B^,,  the  ray  origin  p  lies  in  the  unbounded  face  of  A{B^,). 
Moreover,  observe  that  the  segments  of  B^^  are  contained  in  A^^^ ,  which  imphes  that  p 
hes  in  the  unbounded  face  of  A{B^j ),  for  any  child  of  v.  Therefore  p  Hes  in  the  unbounded 
face  of  ^(Uj<M  ^w  )  (where  M  is  the  number  of  children  of  v).  Hence, 


$(  U    B^,  p)  =  *(JP-„,p). 


uew^ 


The  desired  ^{Q^,p)  is  now  obtained  using  (6.1)  and  Tsi^v). 

6.3      Analysis  of  the  algorithm 

The  correctness  of  the  algorithm  follows  from  the  above  discussion,  so  we  only  have  to 
analyze  the  time  and  space  complexity  of  the  algorithm.  First  consider  the  query  time 
Q{n).  Let  n^  be  the  path  followed  by  the  algorithm  as  it  computes  ^{Q,  p).  We  bound 
the  time  spent  at  each  node  v  G  Tip.  We  spend  O(logr)  time  to  find  the  trieingle  A( 
containing  the  ray  origin  p.  It  follows  from  [EGH*]  (see  also  [Aga],  [Age])  that  ^{C'^^.p) 
can  be  computed  in  0{Jn'^^  log'^'u,,)  ^i™^-  For  each  of  these  triangles  A^,  intersected  by  p, 

we  spend  C>(logn)  time  to  compute  ^{A^,  p)  [CGa];  we  spend  0{Jn^a{ny)\o^  n^)  time 
to  compute  $(,^„,p)  because  Tv  has  at  most  0(n„Q(nt,))  segments  (cf.  Theorem  3.10). 
Therefore,  the  time  spent  at  a  node  v  is  at  most  O(y'n„a(n„)log^  n^  +  \/n^|logn„  + 
5logn„).  Since  s  <  r^  =  0(1)  f  d  n^^  <  n'^,  the  totzd  time  spent  at  the  node  v  is 
0(wnva(n„)  log^  n  +  ^/n^logn).  S   .nming  over  all  nodes  of  lip,  we  obtain 

Q{n)     =      Y.  oUn,a{n,)\og^n  +  ^^ogn]  (6.2) 
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It  is  easily  checked  that,  for  a  node  i'  at  level  i  of  T,  we  have  n^,  <  — .    Moreover  if  :  is 

r' 

the  father  of  v  then,  by  construction,  Aj,  meets  at  most  ^  segments  of  ^,,  which  implies 

K  <^-  Therefore 

Tiy  <  —     and     n'    <  —  (6.3) 

r'  r' 

Thus,  (6.2)  becomes 


Qin)      =     'fol^l^log^n 

1=0 


=     0(\Jna(n)log^  n)  (because  r  >  2) 

Next,  let  us  analyze  the  space  complexity  S{n)  and  the  preprocessing  time  P{n)  of  our 
algorithm.  At  each  node  u  €  T  we  store  the  following  data  structures: 

(i)  Mv-  The  node  v  is  partitioned  into  0(r'^)  triangles  in  0(ni,r  log n^,  log'^~V)  time 
[Agb],  therefore  by  [EGS],  Mv  can  be  preprocessed,  in  time  0{r^  logr),  into  a  data 
structure  of  size  0{r^)  for  point  location  queries.  Since  r  is  chosen  to  be  constant, 
the  time  bound  is  just  Oin^  logn^)  and  the  space  required  is  0{ny). 

(ii)  Ti(£'^,):  It  follows  from  the  result  of  Edelsbrunner  et  aJ.  [EGH*]  (see  aiso  [Age])  that 
Ti(£(,)  requires  0{n'^,  log^  n[)  space,  and  can  be  constructed  in  0{n'l'^  log""  n[)  time. 

(iii)  T2(.4v):  Since  Tiy  has  at  most  0{{nv  +  n[)a{ny  +  n[))  edges,  T2(.4v)  requires  only 
0{{nv+n[)a{ny-irn'^)))  space  (cf.  [CGa]),  and  can  be  constructed  in  0({nv  +  n[)a(nv  + 
Olog(nv  +  n'J)  time. 

(iv)  T3(jFi,):  ^v  has  0(nva(n„))  edges,  therefore  by  Theorem  3.10,  Tji^^)  requires  only 
0(nt,Q:(nv)log"'nv)  space,  and  can  be  constructed  in  0{nl^^a^^^{nv)\og^  n^)  time 
(which  subsumes  the  time  0(nva(n„)  log^  n„)  needed  to  compute  Tv,  as  in  [EGSh]). 

Thus,  the  space  used  at  v  is  0{nva{ny)  log''  n^  +  n[  log^  n^).  Summing  over  all  nodes  of  T, 
we  get 

•^(n)   =    X^(n„a(n  J  log^n^  +  <log^O. 

Observe  that  each  triangle  of  M^  meets  O(^)  segments  of  Q^,  therefore  if  W^  is  the  set 
of  children  of  v,  then 

E   <   =   0{n.r),  (6.4) 

which  implies  that 

•^(")      =      E  0(n„a(n„)  log^  n^  +  n^r  log^  n^) 

=     E  0(n„a(n„)log^nJ 
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Since  each  endpoint  of  a  segment  e  G  ^  falls  inside  one  triangle  A^,  for  each  level  of  T,  e 
is  associated  with  at  most  two  nodes  at  each  level.  Let  l{v)  denote  the  level  of  the  node  v 
in  T.  Then  for  every  i  <  log  n,  we  have 

X]   "v    <   2n.  (6.5) 

l(v)=, 

Hence, 

logn 

S{n)      =      J3  0(nQ(n)log^n)       =       C)(na(n)  log"  n) 
t=i 

Finally,  we  bound  the  preprocessing  time  P{n)  of  our  algorithm  The  above  discussion 
implies  that  the  total  time  spent  in  preprocessing  is  at  most 

P{n)      =      ^0(n^,/V/^(nJlog-n  +  n'//V/Vjlog-n) 

=     ^0(n^/V/2(n,)log-n)  (6.6) 

Since  u)  is  a  constant  <  4.33  and  we  do  not  know  its  exact  value,  we  can  replace  the  term 
a^/^(n„)log'^n  in  (6.6)  with  log"' n  without  any  real  loss  of  any  information.  Therefore 

^(")      =      E    E   C'(n^^^log-n) 
1=1  /(i')=i 

=      EO     r'(^         log-'n 


.=  1  V         ^r'J 


/  log"  1  V 

=  o(„3/=.og-„i:^) 

<     0{n^^^  log"  n)         (because  r  >  2) 
Hence,  we  caxi  conclude  that 

Theorem  6.2  Given  a  collection  Q  of  n  (possibly  intersecting)  segments,  we  can  prepro- 
cess  Q,  in  time  0(n^/^  log'"  n),  into  a  data  structure  of  size  0(na(n)  log'*  n)  so  that,  for 
any  query  ray  p,  we  can  compute  ^{Q,p)  m  0{^Jna{n)\o^  n)  time. 

D 

Remark  6.3:  If  Q  contains  unbounded  segments,  then  the  trizingle  At,  associated  with 
the  root  u  of  T  should  be  a  triangle  that  contains  all  intersection  points  and  all  bounded 
segments  of  Q.  Such  a  A^  can  be  easily  computed  in  C>(nlogn)  time.  Now  for  each 
segment  e   G  ^,  we  compute  e'  =  e  D  A„  and  apply  our  algorithm  to  the  new  set  of 
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seo-ments.  The  portions  of  the  segments  lying  in  the  exterior  of  A^  do  not  intersect  each 
other,  and  are  ordered  in  the  non-decreasing  order  of  their  slopes  along  5Au  in  counter- 
clockwise direction.  Therefore  if  a  query  ray  does  not  hit  a  segment  of  Q  inside  A^.  we  can 
determine,  in  additional  0(log7i)  time,  the  first  segment  hit  by  the  ray  outside  A^,  which 
shows  that  our  algorithm  works  for  unbounded  segments  as  well. 

6.4      Tradeoff  between  space  and  query  time 

In  this  subsection  we  establish  a  tradeoff  between  space  and  query  time  for  ray  shooting 
in  general  arrangements  of  segments.  As  in  Section  4  we  first  give  a  very  simple  algorithm 
that  preprocesses  Q.  in  time  0(n^a^(n)  logn),  into  a  data  structure  of  size  0{n^a^{n))  so 
that,  given  a  query  ray  p,  ^(Q,p)  can  be  computed  in  O(logn)  time. 

Compute  the  arrangement  A{Q)  in  time  0(n^  log  n)  using  the  line  sweep  method  [PS] 
(or  in  time  0{n^)  using  a  more  involved  algorithm  [EOS]),  and  preprocess  A{Q)  for  point 
location  queries  [EGS],  [ST].  Since  the  edges  of  A{Q)  axe  non-intersecting,  we  can  prepro- 
cess each  face  of  A{Q)  into  a  data  structure  T/  for  logarithmic-time  ray  shooting  queries, 
using  0(|n/|^)  space,  where  n/  is  the  number  of  edges  bounding  /,  as  described  in  Section  4. 

To  compute  $(^,p),  for  a  query  ray  p,  first  locate  the  face  /  of  A{Q)  containing  the 
ray  origin  p.  Obviously  ^{Q,  p)  lies  on  the  boundary  of  /  and  therefore  ^{Q,  p)  =  <J*(5/,  p) 
can  be  computed  in  O(logn)  time,  using  T/.  Thus,  the  overall  query  time  is  O(logn). 

As  for  the  storage,  A{Q)  can  be  preprocessed  for  point  location  queries  using  0['n}) 
space  (cf.  [EGS],  [ST]).  The  total  space  required  to  store  all  T/  is  0(^1/  n^})-  Theorem  4.3 
impUes  that  the  preprocessing  time  is  0(^y  n^logn).  It  has  been  shown  in  [EGP*]  that 

^    n]  =  0{n'a\n)). 

Hence,  we  have 

Theorem  6.4  Given  a  collection  Q  of  n  segments  m  the  plane,  we  can  preprocess  Q ,  m 
time  0(n^Q^(n)logn),  into  a  data  structure  of  size  0{n'^a^{n))  so  that,  given  a  ray  p, 
^{G,p)  can  be  computed  in  O(Iogn)  time. 


Next  we  give  an  algorithm  for  the  general  case,  where  n^'*'^  <  rn  <  n^~'',  for  some 
constants  cq,  Ci  >  0.  We  thus  write  m  =  n'',  for  some  l  +  eo^7<2  —  ei.  To  preprocess 
Q  into  a  data  structure  of  size  0{m),  we  proceed  in  the  same  way  as  in  Section  6.1  except 
that  at  each  node  v  we  are  allowed  more  space,  so  we  construct  laxger-size  data  structures 
that  faciUtate  faster  ray  shooting  in  C[,  T„,  etc. 
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Edelsbrunner  et  al.  [EGH*]  (see  also  [Age])  have  shown  that,  given  a  set  £  of  n  Hnes  and 
a  parameter  I  <  di  <  n,  C  can  be  preprocessed,  in  0{n'^^^ s/Bi  log""  n)  time^.  into  a  data 
structure  Ti(£)  of  size  0(nl3i  log'^  n)  so  that,  for  any  query  ray  p,  we  can  compute  $(£,p) 
in  0(,/j-logn)  time.  At  each  node  v  of  level  i,  we  store  Ti(£'y)  with  an  appropriate  value 
of  ;3i  =  /3J  (to  be  specified  later). 

Similarly,  we  have  shown  in  Section  4  that,  given  a  set  £  oin  non-intersecting  segments 
and  a  parameter  /52i  we  can  preprocess  S,  in  time  0[n^'^ s/^log"^  n),  into  a  data  structure 
T3(£')  of  size  0(n/32  log"' n)  so  that,  for  a  query  ray  /?,  we  caji  compute  <^(5,p)  in  time 
0{J^\o^  n).    For  a  node  v  at  level  i,  we  store  T^iTy)  with  aji  appropriate  value  of  ^2 

=  '3\.   (Recall  that  if  Ty  ha^  a  vertex  of  degree  >  3,  then  the  segments  of  Ty  need  to  be 
modified,  as  described  in  Section  4.3.) 

Besides  these  two  types  of  data  structures,  we  also  store  T2(.4v)  of  Section  6.1  at  each 
node  V  oi  T .  For  answering  a  query,  we  use  the  same  procedure,  as  described  in  Section 
6.2. 

Next,  we  analyze  the  complexity  of  this  cdgorithm.  First,  consider  the  space  used  by 
our  algorithm.  Since  \C'^\  =  n^,  \Ty\  =  0{nya{ny))  and  |A^|  <  n'^  +  Uy,  the  space  used  by 
a  node  v  of  T  at  level  i  is 

0{n'J\  log^  n[  +  nya{ny)0'^  log'  n,  +  (n,  +  n[)a{ny  +  n',))  =  0{n'J\  log^  n[  +  n,^^  log'  nj. 

Therefore  the  total  space  used  is 

1=0  /(v)=i 

=      E  0((  E  O^;  V<)  +   E0({  E  ny)|3',a{n)\og'n^ 
t=o       ^  ;{v)=.  ^         1=0       ^  /(v)=.  ^ 

By  (6.4)  and  (6.5),  we  obtain 

Sin)      =      J2  O{0{n\og'n:)  +  xfo(^Jna(n) log' n,) 

1=0  1=0 

Let  n,  (resp.  n',)  denote  the  maximum  value  of  l^^l  (resp.  \Q'J)  for  a  node  v  at  level  i.  By 

(6.3),  we  have  n,,n',  <  ^.  If  we  choose  0{  =       \       and  0^  =  \ — 5 — ,  which  is  easily 

log  n,  a(n)  log  n, 

seen  to  satisfy  0{,  02  >  1,  then  we  obtain 

^^"^      =      ^^(dv"l°g'"0  +  EO    --^^-^-a(n)log'n.) 
1=0      ^^og  ".  ^         ,=0      Va(n)  log'' n,  / 


'Actually  the  preprocessing  time  is  O((n0i  +  n^/^v^log"  n),  but  it  can  be  verified  that  for  our  choice 
of  01  the  first  term  never  dominates,  so  for  simplicity  we  only  write  the  second  term. 
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logn  ^^ 


=    0(1:,,.^).  0(1:;:^ 


logn 


=     0{n'')     because  7  >  1 
=      Oim) 

As  for  the  query  time, 


logn 


n. 


Q(r^)    =    Eo(^/^iogn;)  +  i:o 


logn 


1=0        \  y  "1  /        t=o 

Substituting  the  values  of  l3\  and  d^,  we  obtain 

logn 


n,a(n,)        2 


logn 

Q{n)      =      Y.0 
1=0 


n. 


n,a(n,) 


n;  /         ,=0       V^  "•      /(a(n)log   n,) 


n, 


.VTVlog' 

'°8"         //„xl— r/2  \ 

=      E0((^)  c.(n)log^/^nj 

=     0(n'-^/2a(n)log'/2n) 

Finally,  the  time  spent  in  preprocessing  a  node  u  of  T  at  level  i  is 

Using  the  same  argument  as  for  (6.6),  we  cam  ignore  the  term  a"'/^(n„)  in  the  above  equality. 
Therefore  summing  over  all  nodes,  we  get 

^(")      =      L    E  OK'^'T^log-^  +  n^/^y^log-n,) 

(Substituting  the  values  of  /?J  cind  j3\,  and  using  (6.3)) 
=     0{nyJTp\og^-^l^n) 
=     0{n^\og^-^l^  n) 

Hence,  we  caui  conclude  that 

Theorem  6.5  Given  a  set  Q  of  n  segments  and  a  parameter  n^"*"*"  <  m  <  n^"'' ,  for 
some  constants  eo,ei  >  0,  we  can  preprocess  Q,  xn  0{nyJrn\o^~^^'^  n)  time,  into  a  data 
structure  of  size  0{Tn)   so  that,  for  any  query  ray  p,    we  can   com,pute  ^{Q,p)   in  time 
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n 


Remark  6.6:  The  algorithm  of  [EGH*]  actually  constructs  Ti(£),  in  time 

Oin3i  logn  •  log'^-^  /?:  +  n^'^^^log^  ^), 


using  0(n^i  log  j-)  space,  and  answers  a  query  in  time  0{Jj^  log  -^).  Similarly  the  algo- 
rithm described  in  Section  4  constructs  T3(^),  in  time 

0(n/32  logn  •  log-^  02  +  n^'^Mlog^  ^), 

using  0{n32  log  j-)  space,  and  emswers  a  query  in  0{J^ log  ^)  time.  Using  these  bounds 
in  the  above  analysis,  we  can  improve  the  query  time  Q{n)  to 

0f!^logV^(^)+logn^ 
The  preprocessing  time  is  now  Of  mlog'^n  +  nQ{n)y/m\og'^~^'^{^^^l^)] . 

7     Implicit  Point  Location 

The  planar  point  location  problem  is  a  well  studied  problem  in  computational  geometry 
[Ki],  [EGS],  [ST].  In  this  problem  one  is  to  preprocess  a  given  planax  subdivision  so  that, 
for  a  query  point,  we  can  quickly  determine  the  face  of  the  subdivision  conteuning  it. 
Guibas  et  al.  [GOS]  have  considered  a  generalization  of  this  problem,  in  which  the  map 
is  defined  as  the  arrangement  of  n  possibly  intersecting  polygonal  objects  of  some  simple 
shape,  and  the  goal  is  to  compute,  for  a  query  point  p,  certain  information  related  to  its 
position  within  the  arrangement  of  the  objects;  for  example  to  determine  whether  p  lies 
in  the  union  of  the  objects.  For  simplicity  we  break  the  given  objects  into  a  collection  of 
segments,  and  consider  the  following  formal  statement  of  the  problem: 

We  are  given  a  collection  ^  =  {ci,  . . .  ,  e„}  of  n  segments,  and  with  each  seg- 
ment e  we  associate  a  function  tpe  defined  on  the  entire  plane,  which  assumes 
values  in  some  associative  and  commutative  semigroup  S  (denote  its  operation 
by  +).  Define  ^(i)  =  Yle^s'^Pei^)-  We  want  to  preprocess  Q  so  that,  for  any 
query  point  p,  we  can  quickly  compute  ^(p). 

We  assume  that  rpe  and  ^  satisfy  the  following  conditions: 

(i)  For  amy  given  point  i,  rpe{x)  can  be  computed  in  0(1)  time. 
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(ii)  Any  two  values  in  S  can  be  added  in  0(1)  time. 

(iii)  Given  a  set  ^  of  n  segments  in  the  plane,  we  can  preprocess  it  in  time  0(n  log''  n). 
for  some  constant  k  >  0,  into  a  linear-size  data  structure  'DiQ)  so  that,  given  a  point 
X  lying  either  above  all  the  lines  containing  the  segments  of  Q.  or  below  all  these 
lines,  ^(x)  can  be  calculated  in  O(logn)  time. 

It  is  shown  in  [GOS]  that  many  natural  problems  including  the  problem  of  determining 
whether  p  lies  in  the  union  of  the  given  objects,  or  of  counting  how  majiy  objects  contain 
p,  fall  into  this  scheme.  See  also  the  following  section  for  details. 

The  goal  is  to  come  up  with  an  algorithm  that  uses  O(nlog  *^'  n)  space  and  computes 
^(p),  for  any  query  point  p,  in  sublinear  time.  Guiba^  et  al.  [GOS]  gave  a  randomized 
algorithm,  with  0(n  log''"'"^  n)  expected  running  time,  to  construct  a  data  structure  of  0(n) 
size  so  that,  for  a  query  point  p,  ^(p)  can  be  computed  in  0{n^^^^^)  time,  for  any  6  >  0. 
In  this  section  we  present  an  algorithm  that  improves  the  query  time  to  Oi^/nlog^  n),  and 
makes  the  preprocessing  deterministic  (albeit  no  longer  close  to  linear). 

Let  C  denote  the  set  of  lines  containing  the  segments  of  Q.  Dualize  the  lines  of  C  to 
obtain  a  set  £'  of  n  points.  Let  C  =  {Ci ,  ...  ,Ck}  denote  a  family  of  fc  =  0( log  n )  spanning 
paths  on  £',  with  o-(C)  =  0{y/n).  We  show  how  to  preprocess  a  single  path  C  G  C. 

First,  construct  a  bineo-y  tree  B  =  B{C)  as  in  Section  2.  With  each  node  v  of  B  we 
eissociate  a  set  Q^  of  segments  e  6  ^  such  that  the  dueJ  of  the  line  containing  e  belongs  to 
Sv  (as  defined  in  Section  2).  At  each  node  t;,  we  store  V[Qv)  so  that,  for  any  query  point 
p  lying  either  above  aJI  the  Unes  containing  the  segments  of  Q^  or  below  all  of  them,  ^v(p) 
=  ^egc„  0e(p)  can  be  computed  in  O(logn)  time. 

For  a  given  query  point  p,  we  compute  ^(p)  as  follows.  Let  p*  denote  the  duaJ  of  p. 
Obviously 

Therefore,  it  suffices  to  show  how  to  compute  ^v{p),  for  a  node  v  £  Vb{p').  Observe  that 
for  any  v  G  Vb(p*),  p*  lies  either  above  all  the  points  of  Sv,  or  below  all  of  them,  say  below. 
Since  duality  preserves  the  above-below  relationship,  p  lies  below  all  the  lines  containing 
the  segments  of  Q^-  Therefore,  ^t,(p)  can  be  easily  computed  in  (9(log  n)  time  using  'D{Qv)- 
Next,  let  us  aneilyze  the  complexity  of  our  adgorithm.  First  consider  the  time  spent  in 
answering  a  query.  By  Theorem  2.1,  we  can  determine,  in  O(logn)  time,  a  path  C  e  C 
that  intersects  p*  in  at  most  0{y/n)  edges,  cind  it  follows  from  the  discussion  in  Section  2 
that  Vb{p'),  for  a  given  hne  p',  can  be  computed  in  0{y/nlogn)  time.  By  property  (iii), 
for  each  v  G  Vb{p'),  ^v(p)  can  be  calculated  in  C>(logn)  time.  Thus  the  total  time  spent 
is  0{y/n\og^n).  As  for  the  space  complexity,  V{Qy)  requires  0(|^v|)  space.  Since  the 
segments  associated  with  the  nodes  of  B  at  the  same  level  are  pairwise  disjoint,  the  total 
space  required  to  store  B  is  C>(n  log n).  Finally,  the  preprocessing  time  is  bounded  by  the 
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time  spent  in  computing  C  plus  the  time  spent  in  preprocessing  Qy  for  all  v  E  B.  Hence, 
the  total  preprocessing  time  is  Oin^^"^  log^  n  +  n  log''''"^  n )  =  0(n-^/^  log"' n). 

Therefore,  we  can  conclude 

Theorem  7.1  Given  a  collection  Q  of  n  segments,  and  function  t/v  associated  with  each 
segment  satisfying  properties  (i)-(ui).  one  can  preprocess  Q,  m  0(n''/^  log""  n)  time,  into  a 
data  structure  of  size  O(nlog^n)  so  that,  for  any  query  point  p,  >I'(p)  can  be  computed  m 
0(>/nlog*  n)  time. 


Remark  7.2: 

(i)  As  in  Section  3,  we  can  reduce  the  space  complexity  to  0(n  log  n)  by  maintaining  a 
single  tree  structure  instead  of  a  feimily  of  O(logn)  trees.  Also,  if  we  allow  random- 
ization, then  the  (expected)  preprocessing  time  is  0(n''''"' log^  n),  but  the  query  time 
increases  by  a  factor  of  log  n. 

(ii)  In  some  applications,  where  calculation  of  ^(i)  in  (iii)  above  is  accomplished  by  a 
binary  search,  it  is  possible  to  reduce  the  query  time  to  0(-yn  log  n),  using  fractional 
cascading. 

(iii)  As  in  the  case  of  the  ray  shooting  problem,  the  query  time  cam  be  improved  by  allong 
more  storage.  Instead  of  describing  the  tradeoff  for  the  general  case,  we  will  describe 
it  in  the  next  section  for  a  specific  example. 

(iv)  In  a  companion  paper  [Age]  we  solve  the  batched  version  of  this  problem,  where  ail 
the  query  points  p  are  given  in  advance.  We  present  there  a  solution  that  runs  in 
time 

0{m^^\^^^  log2/3  n  log'/^"  -^  +  n  log*  n  log  -^  +  m  log  n), 

where  m  is  the  number  of  given  query  points. 

8      Other  Applications 

In  this  section  we  consider  other  applications  of  our  technique.  All  these  problems  were 
studied  in  [GOS],  who  obtained  algorithms  with  0(n^/^''"*)  query  time,  for  any  6  >  0.  We 
show  that  using  our  approach  the  query  time  can  be  reduced  to  roughly  y/n. 

8.1      Polygon  containment  problem  —  preprocessing  version 

First  consider  the  following  problem: 
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Given  a  set  Tofn  (possibly  intersecting)  triasigles.  we  wajit  to  preprocess  T  so 
that,  given  a  query  point  p.  we  can  quickly  count  the  number  of  triangles  in  T 
containing  p  (or  just  determine  whether  p  lies  in  the  union  of  these  triangles). 


Figure  15:  Polygon  containment  problem 


We  first  present  an  cdgorithm  that  uses  roughly  linear  space,  and  then  show  that  the 
query  time  caji  be  improved  by  using  more  space.  Our  algorithm  is  based  on  the  following 
observation  of  [GOS].  Let  Q  denote  the  set  of  edges  bounding  the  triangles  in  T  and,  for 
each  e  E  Qy  let  B{e)  denote  the  semi-infinite  trapezoidal  strip  lying  below  e.  Define  a 
function  ip^  in  the  plane  so  that  xpe(p)  =  0  for  a  point  p  outside  Bf.,  and  for  p  E  B^,  U'fip)  = 
1,  if  the  triangle  corresponding  to  e  lies  below  the  line  containing  the  segment  e,  otherwise 
rl)g{p)  =  -1.  It  can  be  checked  that  '^{p),  for  a  point  p,  gives  the  number  of  triangles  of  T 
containing  p.  Moreover,  t/jj  obviously  satisfies  properties  (i)-(ii).  As  to  property  (iii),  if  a 
point  p  lies  above  all  Unes  containing  the  given  edges  then  ^(p)  =  0,  by  definition.  On  the 
other  hand  if  p  hes  below  all  these  lines,  we  do  the  following.  Let  e  denote  the  x-projection 
of  an  edge  e  of  some  triaoigle.  It  is  easily  checked  that 

where  €_,  is  the  non-zero  value  of  V'e,  at  p.  Note  that  the  sum  of  the  right  hand  side  remains 
the  same  between  two  consecutive  endpoints  of  the  projected  segments,  and  the  constant 
values  of  ^  over  these  intervals  can  be  computed,  in  overall  time  0(n  log  n),  by  scanning 
the  projected  segments  from  left  to  right.  Hence,  we  can  preprocess  T,  in  time  0{n\og  n), 
into  a  data  structure  T>  so  that,  for  a  point  p  lying  below  all  lines  of  £,  'I'(p)  can  be 
computed  in  O(logn)  time. 

Thus,  the  observation  of  [GOS]  and  Theorem  7.1  imply  that  by  preprocessing  Q^  into 
the  above  data  structure  V^,  for  each  node  v  oi  B,  the  number  of  triangles  in  T  containing 
a  query  point  p  can  be  counted  in  0(v/nlog^n)  time.   But  observe  that  each  of  the  data 
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structures  Vy  is  a  sorted  list,  and  at  each  node  v  we  do  a  binary  search  in  Vy  to  compute 
<i>y.  We  can  therefore  apply  fractional  cascading  technique  of  [CGb]  to  the  collection  of 
hsts  Vy  attached  to  the  nodes  v  of  5.  This  will  allow  us  to  search  through  the  lists  V^  of 
all  nodes  v  €  Vsi^^)  in  overall  time  O(logn  +  1^0(^)1)  =  0{y/n\ogn).  Hence,  we  have 

Theorem  8.1  Given  a  set  T  of  n  triangles  in  the  plane,  we  can  preprocess  T,  m  time 
0{n'^^^  ■  log'^'n),  into  a  data  structure  of  size  0(n  log  n)  so  that,  given  a  query  point  p,  we 
can  ietermme,  m  time  0{\/n\o^  n),  the  number  of  triangles  m  T  containing  the  point  p. 

D 

We  next  establish  a  tradeoff  between  space  and  query  time  for  the  polygon  containment 
problem.  If  we  allow  0{n^)  space,  then  we  can  construct  the  entire  arrangement  Ti.  of 
UeeC  -^e-  ^^  ^^  easily  seen  that  the  value  of  ^  does  not  change  within  a  face  of  7i.  and 
while  constructing  Ti  we  can  compute  'i>  for  each  of  its  face.  Now  given  a  point  p,  we  can 
compute  ^(p)  in  O(logn)  time  by  locating  p  in  Ti..  Thus  if  we  allow  quadratic  storage, 
the  query  time  can  be  reduced  to  O(log  n).  Next  we  give  an  algorithm  for  the  general  case 
when  n  log^  n  <  m  <  n^. 


Figure  16:  Triajigle  A,  and  segments  of  Q:  SoUd  lines  are  Q,\  dashed  lines  are  Ql 

Let  r  denote  the  set  of  lines  bounding  the  trapezoidaJ  strips  Be,  that  is,  the  lines  con- 
taining the  segments  of  Q  and  the  vertical  lines  passing  through  the  endpoints  of  segments 
in  Q.  Partition  the  plane  into  M  =  0{r'^)  triangles  Ai,  ...  ,Am,  each  meeting  at  most 
^  hnes  of  T  [Agb].  With  each  A,  we  associate  a  set  Q,  consisting  of  aJl  segments  e  e  Q 
such  that  either  e  or  one  of  the  two  downwaxd-directed  vertical  rays  emanating  from  its 
endpoints  intersects  A,  (see  figure  16).  Let  g[  =  Q  -  ^,.  We  can  compute  Q„  for  each  i, 
in  total  time  O(nrlogn).  Since  A,  does  not  intersect  the  boundary  of  B^  for  e  €  Q'„  V'e 
remains  constant  over  A,.  Moreover,  Y.eeC'^'^  f°r  every  A,,  can  be  computed  in  0{nr) 
time,  as  described  in  [Age].  We  preprocess  ^.  into  a  data  structure  of  size  0(Mog^  7), 
using  the  method  just  mentioned.  For  answering  a  query,  we  first  locate  the  triangle  Ak 
containing  the  query  point  p.   Once  we  know  A^,  EeeC  V'e(p)  can  be  determined  in  0(1) 
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time,  and  Hegc*  Ve{p)  can  be  computed  as  described  above.  Since  \Q,\  =  0{^),  the  query 
time  is  ^    ,_  . 

Q{n)    =    O     ,/-log-  +  logn     . 
\  V  r  r  / 

We  need  O(r^)  space  to  store  the  planar  map  formed  by  Aj,  .  .  .  ,  Am  and  0{  ^  log"  ^)  to 
store  the  data  structure  constructed  for  each  Q,.  Therefore,  the  toted  space  used  is 

S{n)      =      0{r^)  +  0(r'  --log' - 

\        r  r 

=     O(nrlog^-). 
r 

If  we  choose  r  =  5 ,  which  is  easily  seen  to  satisfy  1  <  r  <  n,  then  Sin)  =  0(m). 

and  the  query  time  is 


Q{n)      =      O 


^ 


n  ,         "         , 

^ log  — ;=  +  log  n 

m/(n  log^^)      ^^  ^ 


=     0(  ^log^^  +  logn). 


Finally,  the  preprocessing  time  is 


P{n)      =      0{nr  log  n  log""'  r)  +  O  (r^  •  (-]       log"  -  j 


=     Ofm  log- n  +  n^/2     1^^ log"  ^"l 


n 


=     0(m  log"  n  +  n\/m  log"~^ —T=). 


Hence,  we  caji  conclude 

Theorem  8.2  Given  a  collection  T  of  n  (possibly  intersecting)  triangles  in  the  plane,  we 
can  preprocess  T,  in  time  O(mlog"  n  +  n^/m  log"" ^  ^),  into  a  data  structure  of  size  0{m) 
so  that,  for  a  query  point  p,  we  can  count  the  number  of  triangles  of  T  containing  p  in 
timeO{^log'^  +  \ogn). 

a 

Remark  8.3:  The  batched  version  of  this  problem,  when  all  points  axe  given  in  advance, 
can  be  solved,  using  a  different  technique,  in  time  C>(m^/^n^/^  log^^^  n  log"/^ -^  +  (m  + 
n)logn)  [Age]. 
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8.2      Implicit  hidden  surface  removal 

The  next  problem  that  we  consider  is  the  following  version  of  hidden  surface  removal 
problem: 

Given  a  collection  of  opaque  objects  in  3-dimensional  space,  and  a  viewing 
point  a,  we  wish  to  calculate  the  scene  obtained  by  viewing  these  objects  from 
a. 

The  hidden  surface  removal  problem  has  been  extensively  studied  by  many  researchers 
(see  e.g.  [De],  [MK]),  because  of  its  applications  in  graphics  and  other  area^.  For  the  sake 
of  simplicity  let  us  restrict  our  attention  to  polyhedral  objects,  whose  boundary  T  is  a 
collection  {Ai,  ...  ,A„}  of  n  non-intersecting  triangles.  In  the  case  of  tmpUcit  hidden 
surface  removal,  we  do  not  want  to  compute  the  scene  exphcitly,  but  only  to  preprocess 
them  so  as  to  determine  quickly  the  object  seen  at  any  particular  query  pixel  [CS],  [GOS]. 
In  this  subsection,  we  consider  the  following  special  ca^e  of  the  implicit  hidden  surface 
removal  problem.  Let  T  =  {Ai,  ...  ,  A„}  be  a  collection  of  n  non-intersecting  triangles 
such  that  A,  hes  in  the  plcine  z  =  c,,  where  0  <  Ci  <  C2  <  •  •  •  <  c„  are  some  fixed  heights. 
Preprocess  T  so  that,  given  a  query  point  p  on  the  ly-plajie,  one  caji  determine  the  lowest 
tricingle  A,  hit  by  the  upward-directed  vertical  ray  emaciating  from  p. 

[GOS]  have  given  an  algorithm  for  this  problem  that  uses  randomized  processing  and 
has  0(n^''"^'''*)  query  time,  for  any  6  >  0.  Their  algorithm  first  projects  all  tricmgles 
on  the  xy-plane,  and  then  performs  a  binary  search  through  the  sequence  (Aj,  . . .  ,  Aj^) 
of  projected  triangles  to  find  the  first  index  j  such  that  A*  contains  the  query  point  p. 
Each  step  of  the  binary  search  tests  whether  p  lies  in  the  union  of  some  contiguous  block  of 
projected  triangles,  using  the  polygon  containment  algorithm.  Therefore  the  preprocessing 
step  consists  of  constructing  a  binary  tree  2  on  T  whose  leaves  store  the  triangles  of  T  in 
increasing  height,  eind  each  internal  node  w  is  associated  with  a  set  of  triangles  T^,,  stored 
at  the  leaves  of  the  subtree  rooted  at  w.  For  each  node  w  of  Z,  preprocess  T^;  for  the 
polygon  containment  problem,  using  the  algorithm  described  in  Section  8.1.  It  now  follows 
from  the  above  discussion  that  a  query  can  be  answered  by  following  a  path  tt  in  Z  and 
solving  the  polygon  containment  problem  at  each  node  of  n.  Hence  using  Theorem  8.1,  we 
can  conclude 

Theorem  8.4  The  implicit  hidden  surface  removal  problem,  for  an  ordered  collection  of  n 
triangles  m  three-dimensional  space  can  be  solved  m  0{\/n\og^  n)  query  time,  O(nlog~'n) 
space  and  0{n^^^  log"^  n)  preprocessing. 

D 

Remark  8.5: 

(i)  Recently  severed  cdgorithms  for  other  variants  of  the  implicit  hidden  surface  removal 
problem  have  been  developed;  see  [SML],  [Be]. 
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(ii)  As  in  the  case  of  the  polygon  containment  problem,  the  query  time  can  be  improved 
if  we  allow  more  space.  In  particular,  if  we  allow  0(m)  space,  where  n  <  m  <  n'^. 
then  Q{n)  =  0{-^\og^ -^  +  log^  n)  and  P(n)  =  0{m\og^  n  +  ny^log^''  -^). 

(iii)  We  can  easily  modify  our  algorithm  without  affecting  its  time  complexity  so  that  the 
query  point  p  lies  anywhere  in  R"^,  rather  than  lying  on  the  xy-plane.  We  leave  it 
for  the  reader  to  fill  in  the  details. 

8.3      Polygon  placement  problem 

Finally  consider  the  following  problem: 

Let  P  be  a  k-gon  (not  necessaxily  simple)  and  let  A  =  {Ai,  ...  .An}  be  a 
set  of  n  (possibly  intersecting)  triemgles.  Preprocess  them  so  that,  given  a 
(translated)  placement  of  P,  one  can  quickly  determine  whether  P  intersects 
any  of  the  obstacles  at  that  placement. 

Such  a  situation  arises  in  several  applications  [Chb].  A  special  case,  in  which  P  is 
convex  and  the  triangles  cire  non-intersecting,  ha^  been  widely  studied  (see  e.g.  [BZ],  [CD], 
[Fo],  [LS]).  But  the  best  known  solution  for  the  general  case  is  by  [COS],  who  have  given 
Ein  algorithm  with  randomized  preprocessing  and  0((A:n)^/^"'"'*)  query  time,  for  any  6  >  0, 
by  reducing  this  problem  to  the  polygon  containment  problem.  Using  their  technique,  and 
applying  Theorem  8.1,  we  can  easily  obtain  ' 

Theorem  8.6  We  can  preprocess  A  and  P,  m  0({knf^^  log'^  kn)  time,  into  a  data  struc- 
ture of  size  0(A:nlog  kn)  so  that,  given  a  translated  placement  of  P,  we  can  determine  in 
time  0{vkn\ogkn),  whether  P  collides  with  the  obstacles  at  that  placement. 


Remark  8.7:  The  tradeoff  between  space  and  query  time  described  in  Section  8.1  works 
here  as  well.     Therefore,  if  we  allow  0{m)  space,  where  n   <   m   <   n^,  then  Q(n)  = 

0(-7=log^^=  +  logA:n)  and  P{n)  =  0{m\og'^  kn  +  kny/mlog"^-^  kn). 

9      Conclusions 

In  this  paper  we  presented  efficient  algorithms  for  various  problems  involving  collections 
of  segments  in  the  plane,  using  spanning  trees  with  low  stabbing  number.  Although  these 
algorithms  cire  significantly  faster  than  previously  known  algorithms,  several  interesting 
questions  still  remeiin  unanswered: 
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(i)  The  most  challenging  open  problem  is  to  give  non-trivial  lower  bounds  for  the  ray 
shooting  and  the  implicit  point  location  problems.  Recently  Chazelle  [Che]  showed 
that  if  we  cdlow  only  0{n)  space,  then  a  simplex  range  query  (i.e.  counting  the 
number  of  points  of  a  given  set  contained  in  a  query  triangle)  requires  Q(>/n)  time. 
We  conjecture  that  similar  lower  bounds  hold  for  these  problems  as  well. 

(ii)  A  weakness  of  our  algorithms  is  that,  unlike  the  algorithms  of  [GOS],  the  preprocess- 
ing time  is  not  close  to  hnear.  The  most  expensive  step  in  the  preprocessing  of  our 
algorithms  is  the  construction  of  spanning  trees  with  low  stabbing  number,  which 
runs  in  0(n'^^^  log^  n)  time.  Thus  an  intriguing  open  problem  is  to  come  up  with  a 
close- to- linear  algorithm  for  constructing  such  spanning  trees.  In  fact,  a  determinis- 
tic algorithm  with  0{n*^^)  time  complexity  will  also  be  quite  interesting  (recall  that 
Matousek's  randomized  algorithm  runs  in  roughly  n'*/'^  expected  time). 

(iii)  We  have  shown  in  [Age]  that  in  the  batched  version  of  the  implicit  point  location 
problem  (i.e.  when  aJl  query  points  are  given  in  advance)  the  average  cost  of  comput- 
ing ^(p),  over  aJl  given  query  points  P,  caji  be  made  roughly  n^^^  per  point  while  still 
using  close-to-lineeir  space.  It  is  not  known  whether  we  can  obtain  similar  bounds 
for  the  batched  version  of  the  ray  shooting  problem.  A  more  chcdlenging  problem 
is:  Can  we  answer  a  ray  shooting  query  in  o{y/n)  time  without  increasing  storage 
significantly,  if  we  know  the  lines  containing  the  query  rays  (but  not  their  origins)  in 
advance? 

(iv)  The  algorithm  described  in  Section  5  for  reporting  intersections  between  a  given  set 
of  segments  a  query  ray  works  only  for  a  collection  of  non-intersecting  segments.  Can 
we  generalize  it  for  an  arbitrary  collection  of  segments? 

(v)  Micha  Shaxir  has  asked  the  following  question,  which  is  a  generalization  of  the  poly- 
gon contciinment  problem:  Given  a  set  T  of  triangles,  preprocess  them  so  that,  for  a 
query  segment  e,  one  can  quickly  determine  if  e  is  contained  in  the  union  of  triangles 
of  T.  It  will  be  interesting  to  come  up  with  an  efficient  algorithm  using  spanning 
trees  of  low  stabbing  number. 

(vi)  Finally,  there  remeiins  the  task  of  looking  for  other  interesting  problems  that  can  be 
solved  eflBciently  using  the  spainning  trees  of  low  stabbing  number. 
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